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Quasiequilibrium sequences ol binary neutron stars are numerically calculated in the framework of 
the Isenberg- Wilson-Mathews (IWM) approximation of general relativity. The results are presented 
. . . for both rotation states of synchronized spins and irrotational motion, the latter being considered as 

CO I the realistic one for binary neutron stars just prior to the merger. We assume a polytropic equation 

, of state and compute several evolutionary sequences of binary systems composed of different-mass 

»' ' stars as well as identical-mass stars with adiabatic indices 7 = 2.5, 2.25, 2, and 1.8. From our results, 

we propose as a conjecture that if the turning point of binding energy (and total angular momentum) 
locating the innermost stable circular orbit (ISCO) is found in Newtonian gravity for some value of 
the adiabatic index 70, that of the ADM mass (and total angular momentum) should exist in the 
IWM approximation of general relativity for the same value of the adiabatic index. 
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On ' The final stage of coalescing binary neutron stars is one of the most promising sources of gravitational waves for 
ground-based laser interferometers such as GEO600, LIGO, TAMA300, and VIRGO jlj. It is also considered as one 
of the candidates for short-duration gamma-ray burst sources ^2j. With accurate templates of gravitational waves 
— ..^ , from coalescing binary neutron stars, it may be possible to extract informations from signals observed by the above 
^ ■ interferometers iJ, through the matched filtering technique 4]. Therefore, it is essential to theoretically predict the 
I \ wave form of gravitational radiation as well as to understand the details of the physics of binary neutron stars. With 
^ • such a motivation, many efforts have been invested in this topic (some recent reviews can be found in 0,0]) When 
Pt'' summarizing these researches, it is convenient to separate the evolution of binary neutron stars into three phases. 
^ The first one is the pointlike inspiral in which the separation between the two stars is much larger than the radius of 
one star, and the binary system evolves adiabatically. In this stage, the post-Newtonian expansion works excellently, 
and accurate gravitational wave forms have been calculated 0, H 13 ■ The second phase is the intermediate one or 
hydrodynamical inspiral in which the separation between the two components is as small as a few times of the radius 
of one star. This phase is also inspiraling, the time scale of the orbital shrink being still larger than the orbital period, 
hence the qualifier of quasiequilibrium given to it. What differs this phase from the pointlike inspiral is that the tidal 
deformation of the stars is no longer negligible, making necessary a hydrodynamical treatment, in addition to general 
relativity. The final phase is the merger, in which the two stars coalesces dynamically. The first successful merger 
computations have been obtained by Shibata et al. four years ago [^ llCtlTH . and various groups are currently working 
on this subject d El Q El . 

In the present paper, we focus on the intermediate phase (hydrodynamical inspiral). This phase is interesting 
because one might get informations on the equation of state of nuclear matter [l^ . Furthermore, the most realistic 
initial data for computing the merger are given at the end point of sequences in the hydrodynamical inspiral ^ILj ^■ 



' keisuke@provcncc.c.u-tokyo.ac.jp| 
teric.Gourgoulhon@obspm.fr 

Miller and Sue n pointe d out that the circular orbit condition as an initial data is not a good approximation for corotating neutron star 
binary systems llTlllSl . However, as we describe b elow, th e realistic rotation stat e of binary neutron stars prior to merger is irrotational 
flow, and such initial data have been used in Refs. llOtllll . Furthermore, in Ref. Illl . initial data are given at some separations slightly 
larger than the innermost stable circular orbit. We hope that they do not deviate too much from the realistic inspiral for irrotational 
binary neutron stars. 
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Various approaches to compute binary neutron star models in the hydrodynamical inspiral phase have been consid- 
ered: fi ) (s emi-)analytical 0,11^ and (ii) numerical |^ [2^ 0, in Newtonian gravity, (in) (semi-)analytical 
[2^ l27i I2M I29II and fiy ) numerical in the post-Newtonian approximation, and (v) numerical in the framework of 
general relativity |3^ [s^, 0, l35ll36l IstL |^, 's^, . Binary neutron stars gradually decrease their orbital separa- 
tions as a result of the emission of gravitational radiation. However, it is still impossible to integrate the Einstein's 
equations for the thousands of orbits of the hydrodynamical inspiral. Fortunately the time scale of the orbit shrink 
being much longer than the orbital period, the hydrodynamical inspiral phase can be approximated by a sequence 
of steady-state (quasiequilibrium) configurations. In order to combine general relativity and quasiequilibrium, we 
adopt the conformal-flatness condition for the spatial part of the metric — the so-called Isenberg- Wilson-Mathews 
(IWM) approximation of general relativity [4ll IS fi^ (see for a discussion) . This treatment postpones one of 
the goals of the study of the hydrodynamical inspiral, i.e. the theoretical prediction of the wave form of gravitational 
radiation^. Accordingly, we focus on the evolution of various physical parameters during the hydrodynamical inspiral 
and investigate the nature of the innermost orbit (mass-shedding or orbital instability). 

To get a quasiequilibrium confi gur ation of binary neutron stars, we have to specify the rotation state of the system. 
About ten years ago, Kochaneck Hr\ and Bildsten and Cutler H§\ concluded that the irrotational flow is much more 
realistic than the synchronized rotation, because the shear viscosity of nuclear matter is far too low to synchronize 
the spins of the stars with the orbital motion by the dynamical merging. The formulation for solving quasiequilibrium 
configurations of irrotational binary neutron stars in general relativity has been proposed by several authors [49l Isol 
I51I l52l|. In the present work, we have calculated not only irrotational configurations but also synchronized ones in 
order to exhibit the differences between these two extreme states. Beside the type of rotation, some equation of state 
for neutron star matter must be specified to get a binary model. Although some realistic equation of state arising 
from nuclear physics (see e.g. js^) has to be used, we adopt here the polytropic one for simplicity. Nevertheless we 
vary the adiabatic index in the range 7 = 1.8 to 3 to cover the large range of nuclear equations of state published in 
the literature. Quasiequilibrium sequences of binary neutron stars based on a selection of realistic equations of state 
will be presented in a future work. 

Under the above assumptions, we have computed quasiequilibrium sequences of binary systems composed of neutron 
stars with different masses. In the three known binary neutron star systems (observed as binary pulsars with a large 
mass companion), the two components have almost the same mass: the relative difference between the masses of the 
two stars are 4 % for PSR B1913-hl6 [Hil, less than 0.1 % for PSR B1534-hl2 [Hill, and about 1 % for PSR 
B2127-I-11C [s^ l. Until now, all studies on quasiequilibrium sequences of binary neutron stars in the framework of 
general relativity have dealt with only identical-mass binary systems, except for our previous work |35| (hereafter 
Paper III). In that article we have also presented quasiequilibrium sequences of binary neutron stars with a polytropic 
equation of state. However, we limited ourselves to the adiabatic index 7 = 2, because it was the first step to 
relativistic different-mass binary systems. In the present article, we enlarge the range of adiabatic index, investigating 
softer equations of state (7 = 1.8) and stiffer ones (7 = 2.25, 2.5 and 3). 

The plan of the article is as follows: Section ^1 provides a brief summary of the assumptions, basic equations, and 
numerical methods which we adopt. Some code tests not shown in previous articles of this series are given in Sec. lIIII 
and the numerical results are presented in Sec. II VI Sectionals devoted to the discussions and we conclude the paper 
with a summary in Sec. I VII Through the paper, we use units in which c = G — 1, where c is the speed of light and 
G the gravitational constant. 

II. METHOD 

The details about our method and basic equations have been already given in previous papers of this series (Refs. 
[33| and hereafter Paper I and II respectively). We refer the interested reader to those articles and simply 
summarize here the basic assumptions, the basic equations and our numerical method. 

A. Assumptions 

Firstly, we assume the binary system to be in a quasiequilibrium state. This assumption comes from the fact that 
the time scale ^gw of the orbital shrink driven by the emission of gravitational radiation is longer than the orbital 



^ see however Refs. I45l l46l for the computation of gravitational waves from IWM configurations. 
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period Poib- Indeed the following relation holds: 



orb 



6Mtot/ V V 



where d, Mtot, and n denote respectively the separation between two stars, the total mass, and the reduced mass. 
Hence for an equal-mass system, tew > ^orb for d > 6Mtot- It is well known that apart from shrinking the orbits, 
the reaction to gravitational radiation circularizes them. This implies a continuous spacetime symmetry, called helical 
symmetry [l^ . l49l | and represented by the Killing vector: 

where ft is the orbital angular velocity and d/dt and d/d(p are the natural frame vectors associated with the time 
coordinate t and the azimuthal coordinate ip of an asymptotic inertial observer. 

The second assumption is applied to the stress-energy tensor, which we require to have the perfect fluid form: 

Tf_,„ ^ {e+p)u^Ui, +pg^^, (3) 

where e denotes the fluid proper energy density, p the fluid pressure, the fluid 4-velocity, and g^^ the spacetime 
metric. This constitutes a very good approximation for neutron star matter. 

The third assumption concerns the fluid motion inside each star. We consider both cases of synchronized (also called 
corotating) motion and irrotational flow. Note that while only the latter state can be regarded as realistic _42, i48l | , we 
consider both states in order to exhibit the difference between them and thus investigate the hydrodynamical effects 
on the main properties of the binary system. 

The fourth assumption regards the (zero-temperature) equation of state. We choose a polytropic one: 

p = nrP , (4) 

where n is the fluid baryon number density, and k and 7 are some constants. In the present article, we consider that 
the two stars in the binary system have the same n and 7, even if they have different masses. 

Finally, we assume that the spatial part of the metric is conformall y fia t, which corresponds to the Isenberg- Wilson- 
Mathews (IWM) approximation of general relativity 0,E3jE3 (^see |44| for a discussion). In this approximation, the 
spacetime metric takes the form: 

ds^ = -{N^ - B,B')dt^ - 2Bi dtdx' + A^f,jdx'dx^ , (5) 

N being the lapse function, the shift vector of co-orbiting coordinates, A the conformal factor, and the flat 
spatial metric. The accuracy of this approximation should be checked (see Sec. III. A of Paper I). The comparison 
between the IWM results presented here and the non-conformally flat ones will be performed elsewhere [57l| . 



B. Partial differential equations to be solved 

Under the assumptions listed above, we have derived the equations to be integrated in Paper I. Here we present 
only some outline of them. 

The gravitational field equations have been obtained within the 3+1 decomposition of the Einstein's equation |58j| . 
The constraint equations are treated in the so-called conformal thin-sandwich approach (see for a review) , 

which enables to take into account the helical symmetry of spacetime. The trace of the spatial part of the Einstein 
equation combined with the Hamiltonian constraint results in two equations: 

Ai/ = A-kA^{E + S)+ A^K.jK'^ ~V,vV'(3, (6) 
A/? ^ AtiA^S + ^A^K.jK'^ - i(ViZ^vV + V,/3V'/3), (7) 

where Vi stands for the covariant derivative associated with the flat 3-metric fij and A :— V'Vi for the associated 
Laplacian operator. The quantities v and (3 are defined by := IniV and (3 :— ln{AN), and Kij denotes the extrinsic 
curvature tensor of the t = const hypersurfaces. E and S are respectively the matter energy density and the trace 
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of the stress tensor, both as measured by the observer whose 4-velocity is the unit normal n'^ to the t = const 
hypersurfaces (Eulerian observer): 

E := Tf^^n^'n-', (8) 
We have also to solve the momentum constraint, which writes 



S := A^n^T'^. (9) 



AN' + ^V{WjN^) = -16TrNA'^{E + p)W + 2NA'^K'^Wj{3p - Av), (10) 

where A''* := B' + Q{d/d(py denotes the shift vector of nonrotating coordinates. 

Apart from the gravitational field equations, we have to solve the fluid equations. The equations governing the 
quasiequilibrium state are the Euler equation and the equation of baryon number conservation. Both cases of irrota- 
tional and synchronized motions admit a first integral of the relativistic Euler equation (see Sec. II of Paper I), which 
is written as 

if + z/-lnro+lnr = const., (11) 
where H — Inh is the logarithm of the fluid specific enthalpy: 

/.:=^, (12) 
men 

niB being the mean baryon mass. Fq is the Lorentz factor between the co-orbiting observer (we denote his 4-velocity 
by v^) and the Eulerian observer: 

To := (13) 

and r is the Lorentz factor between the fluid and the co-orbiting observers: 

r := -vf'Uf,. (14) 

For a synchronized motion, the equation of baryon number conservation is trivially satisfied, whereas for an irrotational 
fiow, it is written as 

(HA-^ + VH'^i^ = A^hVJJ'o^iH + CHi^'f^tiH - p) + A'^hU^ViTn], (15) 
where ^ is the velocity potential and ( the thermodynamical coeSicient: 

<:^^. (16) 

a Inn 

Fn denotes the Lorentz factor between the fiuid and the Eulerian observer: 

r„ := -n^M,^, (17) 
and Uq is the orbital 3-velocity with respect to the Eulerian observers: 

f/o = -^. (18) 

Within the IWM approximation, the three Lorentz factors introduced so far have the following expressions: 

Fo = {l-A'f,,U^U^)-'/^ (19) 
F„ = {1 - A^ijWWr'/^ (20) 
F = F„Fo(l - A^fijU'ul), (21) 

where W is the fluid 3-velocity with respect to the Eulerian observer: W = Uq for synchronized binary systems, 
whereas 

for irrotational ones. 
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C. Numerical method 



Before going further, let us recall shortly the numerical method that we employ, including the special technique to 
treat the case of non-integer polytropic indices n = 1/(7 — 1) [23.l60l |. 

Through this series of works, we have developed a numerical code which is based on a multidomain spectral method 
with surface-fitted coordinates 0,|62- This code is built upon the C-I--I- library LORENE j6^. 

One of the merits of considering multiple domains is to allow for boundary conditions at infinity, by compactifying 
the outermost domain. Another merit is that we can separate the region including the matter from the other ones. 
This property results in accurate computations, by setting at the boundary between two domains the discontinuities 
which might occur at the stellar surface. Nevertheless, some infinite gradient in the baryon density appears at the 
surface for stiff equations of state (7 > 2). (Even for 7 < 2, if the polytropic index n is not an integer, higher 
derivatives of den sity diverge at the stellar surface). In order to overcome this, we have introduced a regularization 
technique in Ref. |60j and Paper II and we employ it in the present version of the code. 
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FIG. 1: Relative error in the virial theorem of stationary binary figures in general relativity as a function of the orbital 
separation. Solid, dotted, and dashed curves correspond to relative errors in the three variants of the relativistic virial theorem 
(see text). Left (right) panel is for synchronized (irrotational) binary systems. Both of them are calculated for identical-mass 
binaries with the adiabatic index 7 = 2.5 and the compactness M/ R = 0.16 vs 0.16. 



Relative error in the virial theorem 

M/R=O.I4vsO.I4. Synchronized. 7^2.25 M/R=O.I4vsO.I4. Irrotalional. "|t2.25 

10 I I 



10-^ 




Virial error (mass) 
Virial error (GB) 
Virial error (FUS) 




Virial error (mass) 

Virial error (GB) 

Virial error (FUS) 



1.6 1.8 2 2.2 2.4 2.61.6 1.1 



2.2 2.4 2.6 



FIG. 2: Same as Fig. [Hbut for 7 = 2.25 and M/R = 0.14 vs 0.14. 
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FIG. 3: Same as Fig. [T]but for 7 = 1.8 and M/R = 0.10 vs 0.10. 



III. CODE TESTS 



Numerous tests of tlie numerical code have already been presented in Papers I, II, III and Ref. (25|. Consequently 
we presents here tests only regarding (i) the global numerical error measured by the virial theorem, (ii) the numerical 
error in the first law of binary neutron star thermodynamics, (iii) the comparison with previous numerical solutions, 
and (iv) the convergence to post-Newtonian results at large separation. 



A. Virial theorem 



The virial theorem has proved to be a useful tool to check the global error in numerical solutions for stationary fluid 
systems (see e.g. jHOl). A general relativistic version of the virial theorem has been obtained by Gourgoulhon and 
Bonazzola [g^ for stationary spacetimes. It has been recently extended to binary star spacetimes within the IWM 
approximation by Friedman, Uryu and Shibata The virial relation is equivalent to require 

MaDM = MKomar, (23) 

where Madm is the ADM mass: 

Madm = / ^A^^^dS,, (24) 
and AfKomar IS a Komar-type mass, defined by 

MKomar = ^ / V^NdS,. (25) 

The virial theorem obtained by Friedman, Uryu and Shibata 0| (see also Eq. (5.7) of Ref. ^3|) writes 



VE(FUS) = I \2NA^S + —NA^K^K' + —NA(\/i(3\/' (3 - \7^i^Vi^) 
J I sn ■' 4tt 



d 



X 



0. (26) 



The spacetime generated by a binary system is not stationary, due to gravitational radiation. However, within the IWM approximation, 
the gravitational radiation is neglected in the spacetime global dynamics, so that one is able to recover the virial theorem. 
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By a straightforward manipulation, this integral can be recast in the form of the virial theorem as obtained by 
Gourgoulhon and Bonazzola [6^ : 



VE{GB) 



1 

47r 



A(V,^V P - V - 2Vi/3V v) 



0. 



(27) 



Let us stress that the above identity has been derived by Gourgoulhon and Bonazzola [6J| only for stationary space- 
times and that its validity for IWM spacetimes with helical symmetry has been obtained by Friedman, Uryu and 
Shibata HI. 

As error indicators for our numerical solutions, we have evaluated the quantities 



G 



VE{FUS) 



Madm 
VE{GB) 



Mad 



M 



(28) 
(29) 
(30) 



The results are presented in Figs. ^-[31 as a function of the orbital separation do, defined by the distance between 
the two "center of mass" (see Eq. (128) of Paper I for a precise definition). In each figure, the left (resp. right) 
panel is for synchronized (resp. irrotational) binaries. The vertical axis denotes the relative virial errors defined by 
Eqs. (|2Hl)-(|Sni- One can notice in these panels a sudden change in the slope of the curves around the middle of the 
sequences. This is due to the fact that we use five computational domains for each star and the space around it in the 
case of large separation and four domains for small separation, inducing a change in the accuracy of the computation. 
However, since such discrete changes are only a factor two or three better even if we use five computational domains 
for close separation, we decided to use four domains in order to save computational time. The order of magnitude of 
the virial error is ^ 10^^ through the sequence. This is sufficient to discuss the physics of binary systems accurately. 



B. First law of binary neutron star thermodynamics 

The first law of binary neutron star thermodynamics within the IWM approximation has been derived by Friedman, 
Uryu and Shibata 4^]. It is written as 

SMADM = nSJ, (31) 

where SMadm and dJ are respectively the changes in ADM mass and total angular momentum along a constant 
baryon number sequence, and f2 denotes the orbital angular velocity. Since the relation (|31|l is not enforced in our 
code, we may use it to gauge the numerical error. The relative error \{SM abm — ^SJ)/SMa'dm\ has been evaluated 
for some selected points in the sequences by means of a second order polynomial fitting. It is found that this error is 
smaller than a few 10"'^ for large separations and 1 x 10^^ for medium separations. However, it becomes a few % or 
worse for close configurations and around the turning point {SMadm — 0) of the ADM mass. 



C. Comparison with previous numerical results 

We have compared our results with the irrotational sequences computed by Uryu, Shibata and Eriguchi . Since 
only data at the sequence end points are listed in Ref. [33, we have compared for end point configurations only and 
listed the results in Table In the same table, results of Shibata and Uryu are also given for comparison in the 
case 7 = 2. It is found that the ADM mass agree with each other at the order of 10""^, and that the orbital angular 
velocity as well as the total angular momentum agree at a few %. The orbital angular velocity of our results is a few 
% smaller than that of Uryu, Shibata and Eriguchi because the end point of our sequences is still a few % (in term 
of fi) far from the mass-shedding limit, as discussed below (Sec. lIVX)l . On the other hand, the difference in the total 
angular momentum J seems to arise from some systematic error in Ref. [s^ . Indeed our values for J are smaller 
that of Ref. [s^l, although our binaries have a slightly larger separation. Moreover, the new results given by them in 
Refs. I^ls^ for the case of 7 = 2 agrees with ours at the order of 10"'^, while being slightly smaller than our result 
as expected. 
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D. Comparison with point-mass post- Newtonian computations at large separation 
This last type of test will be presented in Sec. IIV (Jl fsee Figs.l^-I^. 
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FIG. 4: Single static neutron star models in the compactness - gravitational mass plane. The right end of each curve is the 
maximum mass model. The gravitational mass M is rescaled by the polytropic constant k: M = ,.^-1/(2(7-1)) ^ 'pjjg (jots 
correspond to values of compactness considered in the present work. 




FIG. 5: Equatorial to polar ratio of the radial derivative of the enthalpy (cusp indicator) as a function of the orbital separation. 
Left (right) panel is for synchronized (irrotational) binary systems with the adiabatic index 7 — 2.5. Solid, dotted, dashed, 
and long-dashed curves denote the cases of compactness M/R — 0.14 vs 0.14, 0.16 vs 0.16, 0.18 vs 0.18, and 0.20 vs 0.20, 
respectively. 



IV. NUMERICAL RESULTS 



The spirit of the quasiequilibrium approach is to model the (gravitational-radiation driven) evolution of a neutron 
star binary by a sequence of steady-state circular-orbit configurations. Along the sequence the baryon number is 
kept constant, since no matter loss occurs during the evolution of each star, until the end of the inspiral stage. Such 
constant-baryon-number sequences are called evolutionary sequences. 

We have computed numerous evolutionary sequences, for different values of the adiabatic index 7 and various 
compactness. We have used 5 (resp. 4) domains for each star and the space around it in the case of a large 
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Equatorial to polai' ratio of tlie radial derivative of enthalpy 

Synchronized case, 7^2.25 IiTolational case, 7^2.25 
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FIG. 6: Same as Fig. |K| but for 7 — 2.25. Solid, dotted, dashed, and long-dashed curves denote the cases of compactness 
M/R = 0.12 vs 0.12, 0.14 vs 0.14, 0.16 vs 0.16, and 0.18 vs 0.18, respectively. 



Equatorial to polar ratio of the radial derivative of enthalpy 
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FIG. 7: Same as Fig. I^but for 7 = 1.8. Solid, dotted, dashed, and long-dashed curves denote the cases of compactness 
M/R = 0.08 vs 0.08, 0.10 vs 0.10, 0.12 vs 0.12, and 0.14 vs 0.14, respectively. 

(resp. small) separation. In each domain, the number of collocation points of the spectral method is chosen to be 
Nr X Ng X = 25 X 17 X 16 or 33 x 25 x 24, where iV^, Ng, and N^p denote the number of collocation points in 
respectively the radial, polar, and azimuthal directions. 

As in Paper III, we parametrize the difference in baryon mass between the two stars by the values of compactness 
of the stars like [(M/i?)oo,star 1 vs (M/i?)oo,star 2], where (M/i?)oo denotes the compactness of the considered star 
at infinite separation or equivalently that of an isolated spherical star with the same baryon number. Hereafter, we 
abbreviate (M/i?)oo,star 1 vs (M/i?)oo,star 2 as M/R for simphcity. According to recent nuclear equations of state, 
the compactness of a M = 1.4 neutron star is in the range 0.12 < M/R < 0.20, depending on the equation 
of state 5^. We therefore explore this range of compactness for our polytropic equations of state. We display the 
compactness parameters used for each adiabatic index 7 in Fig. 01 For each couple of values (7, M/R), one can see in 
Fig. 21 how far the considered model is from the maximum mass configuration. 

Since we have already shown numerous results for 7 = 2 in Paper III, we focus here on results for 7 = 2.5, 2.25, 
and 1.8. We also consider the case 7 = 3 in Sec. However note that this value is too large for a realistic equation 
of state of neutron star. 
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Equatorial to polar ratio of tlie radial derivative of enthialpy 

Synchronized case, y^2.5 Irrotational case, y^2.5 
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FIG. 8: Same as Fig. I^but for different-mass binary systems compared with an identical-mass binary system. The sohd curve 
denotes the case of identical-mass binary system with compactness M/R = 0.16 vs 0.16. Thin and thick dashed curves are, 
respectively, for the more massive and less massive star of a different-mass binary system with M/R = 0.16 vs 0.18. 



Equatorial to polar ratio of tlie radial derivative of enthalpy 
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FIG. 9: Same as Fig. Hbut for 7 = 2.25 The solid curve denotes the case of identical-mass binary system with compactness 
M/R = 0.14 vs 0.14. Thin and thick dashed curves are, respectively, for the more massive and less massive star of a different- 
mass binary system with M/R = 0.14 vs 0.16. 



A. End points of the sequences 

The start of mass exchange between the two stars marks the end of the existence of quasiequilibrium configurations. 
This mass-shedding corresponds to the appearance of a cusp at the surface of the lightest star (or both stars in the 
equal-mass case). This cusp reveals the vanishing of the gradient of enthalpy in the direction of the companion. Hence 
we consider the equatorial to polar ratio of the radial derivative of the enthalpy 

{dH/dr)^q 

,comp 

{dH/dr)poic 

as a dimensionless cusp indicator: % = 1 at infinite separation, whereas x — when the cusp appears. The variation 
of X along various evolutionary sequences is shown in Figs. |Sl~[3for identical- mass binary systems and in Figs. IS]- 
1101 for different-mass systems compared with identical- mass ones. The horizontal axis presents the orbital coordinate 
separation between the two coordinate centers (maxima of density), c?, normalized by the summation of the coordinate 





(32) 
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FIG. 10: Same as Fig. 11 but for 7 = 1.8 The solid curve denotes the case of identical-mass binary system with compactness 
M /R = 0.10 vs 0.10. Thin and thick dashed curves are, respectively for the more massive and less massive star of a different- 
mass binary system with M/R = 0.10 vs 0.12. 
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FIG. 11: Equatorial to polar ratio of the radial derivative of the enthalpy (cusp indicator) as a function of the orbital separation. 
Left (right) panel is for synchronized (irrotational) binary systems with the compactness M/R — 0.14 vs 0.14. Solid, dotted, 
dashed, and long-dashed curves denote the cases of adiabatic indices 7 — 2.5, 2.25, 2, and 1.8, respectively. 



radii of stars in the direction of their companions (ai and a'l)'^: 



ai 



(33) 



This dimensionless quantity is an indicator of the contact between the two stars: d = 1 corresponds to the two stellar 
surfaces touching each other. Note here that d, ai, and a'l are all coordinate lengths, our coordinates being defined 
by the line element 

From Figs. |31 - [7| one can see that the lines seem to focus toward the point d/{ai + a[) — 1 and = in 

the case of synchronized binary systems. This implies that synchronized binary systems of identical-mass stars in 



In Paper III, Eq. (15) should be replaced by this equation. 
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quasiequilibrium end their sequences by the contact between the two stars. It is worth to point out that the curves for 
all cases of compactness take almost the same tracks, but very slightly deviate from each other: the curves of larger 
compactness are located below those of smaller compactness. On the other hand, in the case of irrotational binary 
systems, the curves seem to reach x = at d > 1. This means that the end point of quasiequilibrium sequences of 
irrotational binary systems is a detached configuration with a cusp point. Let us recall that in Paper II, we have 
found the same result in the Newtonian regime. This tendency is enhanced for the case of larger compactness. Note 
here that in Figs. |31and|Sl the curves with larger compactness in the case of irrotational binary systems turn slightly 
toward the vertical axis around x ~ 0.5. This comes from the numerical error. 

The quasiequilibrium sequences should exist up to x = 0- However, the multidomain spectral method we employ 
cannot handle cusp-like figures because in the surface-fitted coordinate procedure we assume that the stellar surface is 
smooth (see Sec. IV. E of Paper I for a detailed discussion of this point). Indeed our numerical scheme fails to converge 
when X — * 0. Monitoring the convergence by the relative difference 6H between the enthalpy fields of two successive 
computational steps, we hardly reach SH < 10~^ for synchronized binary systems when the orbital separation is 
such that X ^ 0.4. For irrotational configurations, although it is still possible to reach SH ^ 10^^ when the orbital 
separation is such that x ~ 0.4 — 0.5, small but unphysical oscillations of the stellar surface appear. This results in 
the spurious (albeit slight) change of direction near the end of some curves in Figs. [3 and El These tendencies are 
enhanced for stiffer equation of state (larger values of 7) and more compact stars (larger values of M / R) . Therefore, 
we stop the computations around x 0.4 — 0.5. 

We compare the evolution of the cusp indicator x for different-mass binaries with that for identical-mass ones 
in Figs. When extrapolating the curves toward x = 0, the sequences seem to terminate by a cusp at the 

surface of the lightest star, both for synchronized and irrotational systems, while the heaviest star keeps its shape 
close to a spherical one. It is worth to note here that as soon as the two masses differ, even slightly, the sequence of a 
synchronized binary system terminates by a detached mass-shedding limit (Roche limit). Indeed the line for the less 
massive star is always located below that for identical-mass stars which reach x = at d = 1. This implies that the 
line for the less massive star seems to always reach x = at d > 1. 

In Fig. ^] the evolution of x is shown for a fixed value of the compactness parameter, M/R — 0.14 vs 0.14, and 
various values of the adiabatic index 7. It appears clearly that the curves of small 7 are located below those of large 
7. This means that at fixed relative orbital separation d and fixed compactness, a configuration with small 7 is closer 
to the mass-shedding limit than a configuration with large 7. The reason is that stars with small adiabatic index 
have centrally concentrated structures with some extended low density halos. Such low density outer layers result in 
a smaller x- 



Relative change in central energy density 

Synchronized case, y=2.5 Iirotational case. Y^2.5 




FIG. 12: Relative change in central energy density as a function of the orbital separation between the centers of mass of 
each star. Left (right) panel is for synchronized (irrotational) binary systems with the adiabatic index 7 = 2.5. Solid, dotted, 
dashed, and long-dashed curves denote identical-mass binary systems with compactness M/R — 0.14 vs 0.14, 0.16 vs 0.16, 0.18 
vs 0.18, and 0.20 vs 0.20, respectively. 
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FIG. 13: Same as Fig. I12l but for 7 — 2.25. Solid, dotted, dashed, and long-dashed curves denote identical-mass binary systems 
with compactness M/R = 0.12 vs 0.12, 0.14 vs 0.14, 0.16 vs 0.16, and 0.18 vs 0.18, respectively. 
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FIG. 14: Same as Fig. 1121 but for 7 — 1.8 Solid, dotted, dashed, and long-dashed curves denote identical-mass binary systems 
with compactness M/R = 0.08 vs 0.08, 0.10 vs 0.10, 0.12 vs 0.12, and 0.14 vs 0.14, respectively. 



B. Relative change in central density 

Let us now examine the relative change in central energy density along an evolutionary sequence. In Figs. I12l - I14| 
in depicted the variation of the quantity 

Se, (34) 

where Cc and Ccoo denote the central proper energy density of respectively the actual star and an isolated spherical 
star with the same baryon mass (or in other words, that at infinite separation). One can see from these figures that 
for all sequences the central energy density decreases as the orbit shrinks. It is also found that this decrease is larger 
for larger compactness, at fixed adiabatic index and rotation state. Note that the magnitude of the decrease is very 
different between the synchronized case and the irrotational one: ~ 10% for synchronized binaries versus a few % for 
irrotational ones. 

The relative change in central energy density is presented as a function of the cusp indicator x in Figs. 1151 [T7I 
in order to evaluate the total change at the end of the sequence (x = 0). It is found from these figures that the 
central energy density for large compactness decreases more than that for small compactness. Furthermore, we can 
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FIG. 15: Relative change in central energy density as a function of the cusp indicator x- Left (right) panel is for synchronized 
(irrotational) binary systems with the adiabatic index 7 — 2.5. Solid curves with open circle, square, diamond, and triangle 
denote identical-mass binary systems of compactness M/R — 0.14 vs 0.14, 0.16 vs 0.16, 0.18 vs 0.18, and 0.20 vs 0.20, 
respectively. Thick and thin dotted curves are for the less and more massive stars of the binary system with M/R = 0.14 vs 0.16, 
thick and thin dashed curves for those with M/R = 0.16 vs 0.18, and thick and thin long-dashed curves for those with 
M/R = 0.18 vs 0.20. 
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FIG. 16: Same as Fig. I15l but for 7 = 2.25. Solid curves with open circle, square, diamond, and triangle denote identical-mass 
binary systems of compactness M/R = 0.12 vs 0.12, 0.14 vs 0.14, 0.16 vs 0.16, and 0.18 vs 0.18, respectively. Thick and 
thin dotted curves correspond respectively to the less and more massive stars of the binary system with M/R = 0.12 vs 0.14, 
thick and thin dashed curves for those with M/R = 0.14 vs 0.16, and thick and thin long-dashed curves for those with 
M/R = 0.16 vs 0.18. 



see that the decreasing track of central energy density is determined by the compactness of the star itself and not 
by that of its companion, even though the value itself is affected by the companion star. This effect is clearly found 
for synchronized binary systems. For irrotational systems, since the amount of decrease is very small and the curves 
include some numerical errors, it is not clear whether the decrease depends only on the compactness of the star and 
not on that of its companion. However the same tendency could be found, in particular in the case of a soft equation 
of state like 7 = 1.8 (see Fig. I17|) . Note here that the slight changes of direction in the curves around x ^ 0.5 for 
large compactness in the irrotational case in Figs. 1151 and [T^ is due to the numerical error discussed in the previous 
subsection IIV Al 

Although the sequences stop at around x — 0.4, it is possible to predict how much the central energy density 
decreases at the mass-shedding limit by extrapolating the curves in Figs. El^^EI For instance, the central energy 
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FIG. 17: Same as Fig. 1151 but for 7 — 1.8. Solid curves with open circle, square, diamond, and triangle denote identical-mass 
binary systems of compactness M/R = 0.08 vs 0.08, 0.10 vs 0.10, 0.12 vs 0.12, and 0.14 vs 0.14, respectively. Thick and thin 
dotted curves are for the less and more massive stars of the binary system with M/R — 0.08 vs 0.10, thick and thin dashed 
curves for those with M/R = 0.10 vs 0.12, and thick and thin long-dashed curves for those with M/R = 0.12 vs 0.14. 
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FIG. 18: Relative change in central energy density as a function of the cusp indicator x- Left (right) panel is for synchronized 
(irrotational) binary systems with compactness M/R — 0.14 vs 0.14. Solid, dotted, dashed, and long-dashed curves denote the 
cases of 7 = 2.5, 2.25, 2, and 1.8, respectively. 



density of synchronized binary systems with 7 2.25 may decrease by 10 % - 15 % in the range of M / R~ 0.12 vs 0.12 
- 0.18 vs 0.18. 

Figure [THl is shown in order to investigate the behavior of the relative change in central energy density when varying 
the adiabatic index and keeping the compactness fixed. It is found that the rate of decrease is lower for larger adiabatic 
indices. This is because the larger the adiabatic index, the closer the matter to the incompressible state in which the 
central energy density does not change. 



C. ADM mass and total angular momentum 



The variation of the ADM mass A/adm and total angular momentum J along an evolutionary sequence is presented 
in Figs. I19I - I21I These figures show that a turning point (minimum of A/adm or J) appears for synchronized binaries 
with 7 = 2.5 and 2.25. On the other hand, we cannot find it in the synchronized case with 7 ~ 1.8 and in any of 
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FIG. 19: ADM mass and total angular momentum of the binary system as a function of the orbital angular velocity. Top 
(bottom) panels are for identical-mass (different-mass) binary systems of compactness M/R = 0.16 vs 0.16 {M/R — 0.16 vs 0.18) 
with the adiabatic index 7 = 2.5. In each panel, the solid curve denotes the synchronized case, and the dashed one the 
irrotational one. 
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FIG. 20: Same as Fig. 1191 but for 7 — 2.25. Top (bottom) panels are for identical-mass (different-mass) binary systems of 
compactness M/R = 0.14 vs 0.14 {M/R = 0.14 vs 0.16). 



irrotational sequences. Note that the position along the sequence of the turning point in the ADM mass and the total 
angular momentum coincide with each other, as expected from Eq. (|31|l . The turning points are also found in the 
different-mass case for synchronized binary systems with 7 = 2.5 and 2.25. However, their position is closer to the 
end point of the sequence than that in the identical-mass case. Since, as mentioned above, we stop the computation 
slightly before the end point of quasiequilibrium sequences, there remains the possibility that the turning point may 
exist in the range < x < 0.4, even if we do not find it in the synchronized case with 7 = 1.8 and in the irrotational 
case. The complete discussion about the existence of the turning point for various adiabatic indices and compactness 
will be presented in Sec. 

In Figs. 1^21 -ESI we present the ADM mass as a function of the orbital angular velocity, as in Figs.1191- 1211 but with 
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FIG. 21: Same as Fig. CHI but for 7 = 1.8. Top (bottom) panels are for identical-mass (different-mass) binary systems of 
compactness M/R = 0.10 vs 0.10 {M/R = 0.10 vs 0.12). 





-0.002 
-0.004 
1° -0.006 
^ -0.008 
's -0.01 

Q 
_< 

S -0.012 
-0.014 
-0.016 







"1 ' I ' r 




■lfcl.8 M/R=0.08 corot. 
1^1.8 M/R=0.08 irrot. 

- y=2 M/R=0. 12 corot. 

- - Y=2 M/R=0.12 iiTot. 
Y=2.5 M/R=0. 16 corot. 

- - -1^2.5 M/R=0.16 irrot. 
^=2.5 M/R=0.20 corot. 

- - 1^2.5 M/R=0.20 irrot. 

3PN EOB corot. (Damour et al. 2002) 

- - 3PN EOB inot. (Damour et al. 2000) 

3PN non resum. corot. {Blanchet 2002) 

. . ■ . 3PN non resum. irrot. (Blanchet 2002) 
corot. binary BH (Grandclement et al. 2002) 




0.02 



0.04 



0.06 



0.08 



0.1 



0.12 



0.14 



FIG. 22: Relative binding energy along evolutionary sequences of equal-mass binary neutron stars, compared with (i) analytical 
results at the 3rd post-Newtonian order for point-masses by Damour et al. 2000 \6(^ . Damour et al. 2002 and Blanchet 
2002 |6^ . and with (ii) numerical results for corotating binary black holes by Grandclement et al. 2002 [g^. fl is the orbital 
angular velocity and Mo is twice the gravitational mass of a single static neutron star (resp. black hole) of the same baryon 
number (resp. same horizon area) as that defining the considered sequence. 



a scaling defined by twice the gravitational mass of the single static neutron star of the same baryon number as that 
of the sequence, Mq (whereas in Figs. E|- [21 the scaling was set by the polytropic constant k). Note that Mq can 
also be viewed as the ADM mass at infinite separation. This enables us to compare with third order post-Newtonian 
results for point-mass particles obtained in the Effective One Body (EOB) approach by Damour et al. [6^ or 
in the standard non-resummed post-Newtonian framework by Blanchet |68j| . This also enables us to compare with 
corotating binary black holes, according to the numerical results by Grandclement et al. [g^- We note from these 
figures that the irrotational binary neutron star sequences are very close to the 3PN irrotational ones (dashed and 
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FIG. 23: Enlargement of Fig. I22l in the neutron star region, with some additional sequences. 

dotted fine curves). Moreover, all curves converges at large separation (small value of HMq), which can be considered 
as a check of our numerical computations, especially the way fl is determined (cf. Sec. IV. D. 2 of Paper I). 

Another interesting feature shown in Fig. is the good agreement between the irrotational binary neutron star 
sequence with high compactness {M/R = 0.20) and the binary black hole sequence. The latter is made of corotating 
black holes, but for QMq < 0.06, the spin effects are not very important, so that it is meaningful to compare the 
corotating black hole sequence up to fiMp < 0.06 with the irrotational neutron star sequence. Finally, we notice from 
Fig. 1231 that at a given compactness and separation, the configuration with smaller adiabatic index is more bound. 
Conversely, at a given adiabatic index and separation, the configuration with higher compactness is more bound. 
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FIG. 24: Axial ratios a2/a\ (thick curves) or as/ai (thin curves) as a function of the cusp indicator x along sequences with 
the adiabatic index 7 = 2.5. The left (right) panel is for synchronized (irrotational) binaries. Solid, dotted, dashed, and 
long-dashed curves correspond to M/R = 0.14 vs 0.14, 0.16 vs 0.16, 0.18 vs 0.18, and 0.20 vs 0.20, respectively. 
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FIG. 25: Same as Fig. 1241 but for 7 — 2.25. Solid, dotted, dashed, and long-dashed curves correspond to M/R = 0.12 vs 0.12, 
0.14 vs 0.14, 0.16 vs 0.16, and 0.18 vs 0.18, respectively. 
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FIG. 26: Same as Fig. 1241 but for 7 — 1.8. Solid, dotted, dashed, and long-dashed curves correspond to M/R — 0.08 vs 0.08, 
0.10 vs 0.10, 0.12 vs 0.12, and 0.14 vs 0.14, respectively. 



D. Deformation of the stars 



The variation of the stellar shape along an evolutionary sequence is presented in Figs. [^-[2^1 In these figures, we 
depict the (coordinate) axial ratio a2/ai and a^/ai as a function of the cusp indicator x- This allows one to speculate 
about the values of the axial ratios at the end point of the sequence, because the curves are close to straight lines. 
Since x = corresponds to the end point of the sequence (mass-shedding limit), we can evaluate the deformation 
of the stars at the end of the sequence by extrapolating the curves to x — 0. The axial ratios at the end of the 
sequences are thus predicted to be 02/01 ~ 0.7 and 03/01 ^ 0.6 for synchronized binary systems with 7 = 2.5, and 
02/01 ~ 0.67 — 0.72 and 03/01 ~ 0.7 — 0.75 for irrotational ones with 7 = 2.5. It is also possible to perceive that for the 
same value of x the stars with larger compactness are closer to the spherical figure than those of smaller compactness. 

We also present the change of the axial ratios along the sequence in Fig. 1271 by fixing the compactness and varying 
the adiabatic index. It is clearly seen from this figure that that the deformation of the star is larger for larger adiabatic 
index. 



20 



Axial ratios 



M/R=0.14vs0.14, Synchronized M/R=0. 14vs0.14, Irrotational 

1 I ' , ' , ' , ' , ' , ' , ' , ' , ' , . 



0.9 



0.6 







//'/' ' 
/'// 
//■■/ 

//'v 

//// 
'/ 


/AW 
/ /// 

//■V 


/// 


'9/ 


1^2.5 : a>| 


-p2.5 : a,/a| 


7^2.5 : a/:i, 


T=2.5 : a,/aj 


7^2.25 : a^/a. 


1=2.25 : a/a. 


y-2.25 : a>, 


1=2.25 : a/a. 


V=2 : aVa, 


V=2 : a/a. 


Y=2 : a,/a, 


1=2 : a/a. 


1^1.8 : a/a. 


1=1.8 : a,/a, 


7^1.8: a/a. 


1=1.8 : a/a. 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 



X % 



FIG. 27: Axial ratios as a function of the cusp indicator x- Left (right) is for synchronized (irrotational) binary systems with 
compactness M/R — 0.14 vs 0.14. Solid, dotted, dashed, and long-dashed curves correspond to the adiabatic indices 7 = 2.5, 
2.25, 2, and 1.8, respectively. Thick and thin curves are respectively for the axial ratios 02/11 and as/ai. 



E. Contours of baryon density 

Isocontours of the baryon density of various binary systems are shown in Figs. 1281 - 1501 There are two panels for 
different-mass binary systems in each figure. The left panel shows the result for the synchronized case while the right 
one is that for the irrotational case. In each panel, the left-hand side star is the less massive star. It is clearly seen 
that the lighter star is tidally deformed and elongated while the more massive star does not deviate very much from 
the spherical shape. 
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FIG. 28: Isocontours of the baryon density for binary systems with the adiabatic index 7 = 2.5 and the compactness 
M/R = 0.16 vs 0.18. The left panel is for the synchronized case with the orbital separation dc/^poiy ~ 1.356, while the right 
panel is for the irrotational one with da /Rpoiy ~ 1.393. The orbital axis is located on a; = in each panel. The thick solid 
lines denote the surface of the stars. 
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Baryon density (y=0) 

Baryon density (y=0) 




FIG. 29: Isocontours of the baryon density for binary systems with the adiabatic index 7 = 2.25 and the compactness 
M/R = 0.14 vs 0.16. The left panel is for the synchronized case with the orbital separation dG/^?poiy = 1.765, while the right 
panel is for the irrotational one with do /Rpoiy = 1.775. The orbital axis is located on a; = in each panel. 
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FIG. 30: Isocontours of the baryon density for binary systems with the adiabatic index 7 = 1.8 and the compactness 
M/R — 0.10 vs 0.12. The left panel is for the synchronized case with the orbital separation dc/^poiy ~ 4.391, while the right 
panel is for the irrotational one with da/Rpdy — 4.257. The orbital axis is located on a; = in each panel. 



F. Tables of results 



Finally, we summarize our numerical results about constant baryon number sequences in Tables Hll - IXIIll In each 
table, we first give the orbital separation defined as the distance do between the two "center of mass" (see Eq. (128) 
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of Paper I for a precise definitfon), in polytropic units: 

da ■■= (35) 
where -Rpoiy is the length constructed from the the polytropic constant and the adiabatic index 7: 

Rpoiy:=^'/'^'-'^. (36) 

We also list in Tables Hll - IXIIII the relative separation d defined by Eq. 1)33(1 (let us recall that d = 1 indicates the 
contact between the two stars). Next we give the dimensionless orbital angular velocity: 

n:=f]Ki/2(''-i), (37) 

the dimensionless ADM mass: 

M Mk-i/2(^-i), (38) 

the dimensionless baryon mass: 

Mb := Mbk-^^^^^~^\ (39) 
and the dimensionless total angular momentum: 

J -.^ Jk'^/^'-'-^I (40) 

The absolute values of the relative errors in the virial theorem, [^^^^(M)!, \VE{FUS)\, and \VE{GB)\, are defined 
by Eqs. ((28|l - ((3()|l . Sec is the relative change in central energy density defined by Eq. ((34(1 . In the tables, primes 
denote values for the more massive star, and the symbol f denotes the turning point in the ADM mass (and the total 
angular momentum) along the sequence. 



V. DISCUSSION 

As shown by Friedman et al. the minimum (turning point) in the ADM mass in a constant baryon number 

sequence locates the appearance of a secular (resp. dynamical) instability for synchronized (resp. irrotational) binaries, 
thus defining the innermost stable circular orbit (ISCO). The frequency at the ISCO is a potentially observable 
parameter by the gravitational wave detectors, and thus a very interesting quantity. The position of the ISCO for 
identical-mass binary systems is depicted as a function of the compactness parameter in Fig. 1311 In this figure, the 
Newtonian results are shown as M/R = 0. The corresponding values of x for the irrotational binaries are slightly 
different from those in Paper II. This is because we have improved our numerical code to decrease the aliasing error 
which appears in the spectral method [70j . 

Several things can be found and predicted from Fig. |^ Firstly, the turning points are found for synchronized 
binary systems for 7 > 2 and for irrotational ones for 7 > 2.5. This last finding agrees with the results of Uryu et al. 
|37j . The turning point might exist for adiabatic index slightly lower than 2 for synchronized binary systems and 2.5 
for irrotational ones. However, since it is very difficult for our numerical code to calculate quasiequilibrium figures of 
X < 0.4 for mildly relativistic cases and x < 0-5 for highly relativistic ones, we cannot determine the exact value of 
the minimum adiabatic index which has the turning point. Secondly, the quantity x the turning point increases 
proportionally to the compactness of the stars. This implies that the relative orbital separation at the turning point 
increases with the increase of the compactness. We should comment here that the tendency for the irrotational binary 
system with 7 = 3 was formerly given by Uryu et al. |3^ . Thirdly, the inclination of lines for synchronized binary 
systems is larger than that for irrotational ones. Fourthly, we may assert that a turning point is found in Newtonian 
calculations for some value of the adiabatic index 7, it should exist in the general relativistic computations (within 
the IWM approximation) for the same value of 7, because the lines in Fig. |^have a positive slope. It is worth to 
point out that even if the turning point is not found in the Newtonian sequences, the possibility to find it in the 
general relativistic computations does not vanish. In Fig. 1311 some vertical lines which denote the maximum values 
of compactness for stable spherical stars are depicted. They show the limit of the existence of the turning point for 
each adiabatic index. Then, for a given adiabatic index, the turning point could exist from a compactness M/R > 
to the maximum {M / R)^i^x taking into account the fact that the line of the turning point has a positive inclination. 



23 



Turning points 



Maximum values y=l.S I -{=115 y=i 

of compactness (0.|72) , (0-252) (0.pi6) 



for stable spherical stars 




□ 01^2.5 (Synchronized) 

0.2 \- ^ - ^1^2.25 (Synchronized) 
A- —A 7=2 (Synchronized) 
9y=i (Irrotational) 
■ 1^2.5 (hTotational) 



(0.?15) 



(0. 



5 
280) 



0.1 0.2 
Compactness (M/R) 



0.3 



FIG. 31; Positions of the turning points of the ADM mass (and total angular momentum) in the sequence as a function of 
the compactness parameter. Solid curve with filled circle is for the irrotational binary system with 7 = 3. Dotted with open 
square, dashed with open diamond, and long-dashed with open triangle denote the results of synchronized binary systems with 
7 — 2.5, 2.25, and 2, respectively. Vertical solid, dotted, dashed, and long-dashed lines are maximum values of compactness for 
stable spherical stars (at X = 1) for 7 = 3, 2.5, 2.25, 2, and 1.8, respectively. 



VI. SUMMARY 



We have developed a numerical code for computing quasiequilibrium sequences of binary neutron stars through the 
series of works 0, |3, H^- present article, we have presented the results of relativistic computations of 

both synchronized and irrotational rotation states, for polytropic equations of state with several adiabatic indices. 
We have considered not only binary systems composed of identical-mass stars but also different-mass systems. We 
have investigated the behaviors of various physical quantities along constant baryon number sequences (evolutionary 
sequences): the relative change in central energy density, the ADM mass, the total angular momentum, the shape 
of the figures, as well as the location of the end point and that of the turning point in the ADM mass and angular 
momentum. 

As one of interesting results, we propose, as a conjecture, that if the turning points defining the ISCO are found for 
the Newtonian calculations for some adiabatic index 7, they should exist in the general relativistic computations for 
the same value of 7. Of course, this conjecture is derived from the results of non-realistic cases such as synchronized 
binary systems or irrotational ones with 7 = 8. Furthermore, since the inclination of lines in Fig. 1311 is positive, it is 
not possible to exactly predict whether or not a turning point appears in the IWM approximation when there is no 
turning point in the Newtonian calculation. However, if the computations of irrotational binary systems for x < 0-2 
could be obtained in the future with sufficient accuracy in Newtonian gravity, which is easier to be calculated than 
in the relativistic framework, it would be possible to speculate about the turning points in the IWM approximation 
without computations. 

The numerical results presented in this article are freely available from http://www.lorene.obspm.fr/data/ This 
web page provides binary data files containing the whole binary neutron stars configurations as well as some codes to 
read them. 
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TABLE I: Comparison with the results of Uryu, Shibata and Eriguchi and Shibata and Uryu j4Ml65|. 





AI/R 


Mb 


0.5 X M 




J 1 -^^ADM 


da/M 










7 = 2.5 








this work 


0. 


14 


0.08742 


0.07974 


0.01836 


0.9563 


8.815 


Uryu et al. 


0. 


,14 


0.0874 


0.0797 
7 = 2.25 


0.0189 


0.981 


8.56 


this work 


0. 


,12 


0.09316 


0.08639 


0.01414 


1.004 


10.42 


Uryu et al. 


0. 


,12 


0.0932 


0.0864 


0.0145 


1.03 


10.3 


this work 


0. 


,14 


0.1080 


0.09886 


0.01825 


0.9521 


8.812 


Uryu et al. 


0. 


14 


0.108 


0.0988 


0.0189 


0.975 


8.52 


this work 


0. 


,12 


0.1299 


7 = 2 
0.1211 


0.01423 


0.9972 


10.42 


Uryu et al. 


0. 


,12 


0.130 


0.121 


0.0144 


1.02 


10.3 


this work 


0. 


,14 


0.1461 


0.1347 


0.01809 


0.9474 


8.805 


Uryu et al. 


0. 


,14 


0.146 


0.135 


0.0187 


0.971 


8.52 


Shibata & Uryu 


0. 


,14 


0.14614 


0.1347 


0.01885 


0.9419 


8.554 


this work 


0. 


,10 


0.1725 


7 = 1.8 
0.1635 


0.01038 


1.065 


12.88 


Uryu et al. 


0. 


,10 


0.173 


0.164 


0.0106 


1.08 


12.6 


this work 


0. 


,12 


0.1920 


0.1803 


0.01394 


0.9947 


10.49 


Uryu et al. 


0. 


,12 


0.192 


0.180 


0.0143 


1.02 


10.2 


this work 


0. 


,14 


0.2065 


0.1924 


0.01767 


0.9456 


8.871 


Uryu et al. 


0. 


,14 


0.207 


0.192 


0.0185 


0.967 


8.50 



27 



TABLE II: Orbital angular velocity, ADM mass, total angular momentum, virial errors, axial ratios, relative change in central 
energy density, and the cusp indicator x along constant baryon number sequences of identical mass binaries in synchronized 
motion, for various compactness. The symbol f denotes the turning point of the sequences. 
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1 


439 





2077 





1814 


3 


449(-2) 


4.664 


-5) 


8.490 


-5) 


7.876( 


-5) 





8894 





8364 


-5 


156( 


-2) 





6204 


1 


422 


1 


392 





2139 





1814 


3 


450(-2) 


4.221 


-5) 


7.586 


-5) 


6.990( 


-5) 





8780 





8222 


-5 


580( 


-2) 





5849 


1 


367 


1 


300 





2258 





1814 


3 


457(-2) 


3.944 


-5) 


6.632 


-5) 


6.120( 


-5) 





8510 





7903 


-6 


496 ( 


-2) 





5026 


1 


349 


1 


265 





2301 





1814 


3 


462(-2) 


4.404 


-5) 


6.957 


-5) 


6.392( 


-5) 





8390 





7767 


-6 


862( 


-2) 





4663 




















M/R 


= 0.18 vs 0.18, 


Mb 


= 0.1131 vs. 


.1131 




















1 


822 


2 


019 





1582 





2016 


4 


191(-2) 


2.354 


-5) 


4.687 


-5) 


4.417( 


-5) 





9599 





9342 


-2 


481( 


-2) 





8489 


1 


640 


1 


794 





1824 





2016 


4 


133(-2) 


1.872 


-5) 


3.723 


-5) 


3.505( 


-5) 





9434 





9084 


-3 


410( 


-2) 





7885 


1 


531 


1 


651 





2000 





2015 


4 


110(-2) 


3.294 


-5) 


5.649 


-5) 


4.700( 


-5) 





9277 





8853 


-4 


231( 


-2) 





7339 


1 


459 


1 


551 





2135 





2015 


4 


100(-2) 


3.069 


-5) 


5.196 


-5) 


4.248( 


-5) 





9129 





8648 


-4 


958( 


-2) 





6845 


tl 


409 


1 


479 





2235 





2015 


4 


098(-2) 


6.083 


-5) 


1.115 


-4) 


1.021( 


-4) 





8999 





8474 


-5 


568( 


-2) 





6418 


1 


368 


1 


414 





2327 





2015 


4 


100(-2) 


5.309 


-5) 


9.712 


-5) 


8.790( 


-5) 





8858 





8292 


-6 


189( 


-2) 





5967 


1 


313 


1 


322 





2458 





2015 


4 


110(-2) 


4.138 


-5) 


7.013 


-5) 


6.147( 


-5) 





8605 





7986 


-7 


200( 


-2) 





5182 


1 


295 


1 


287 





2505 





2015 


4 


116(-2) 


4.174 


-5) 


6.643 


-5) 


5.739( 


-5) 





8493 





7857 


-7 


606( 


-2) 





4840 




















M/R 


= 0.20 vs 0.20, 


Mb 


= 0.1250 vs. 


.1250 




















1 


822 


2 


140 





1632 





2200 


4 


847(-2) 


3.860 


-5) 


7.632 


-5) 


7.280( 


-5) 





9670 





9442 


-2 


554( 


-2) 





8694 


1 


568 


1 


814 





1994 





2199 


4 


760(-2) 


5.257 


-5) 


9.409 


-5) 


8.084( 


-5) 





9464 





9110 


-3 


973( 


-2) 





7906 


1 


459 


1 


665 





2192 





2198 


4 


737(-2) 


4.707 


-5) 


8.401 


-5) 


7.058( 


-5) 





9310 





8878 


-4 


949 ( 


-2) 





7350 


1 


422 


1 


613 





2266 





2198 


4 


733(-2) 


4.499 


-5) 


8.008 


-5) 


6.666( 


-5) 





9243 





8780 


-5 


357( 


-2) 





7113 


tl 


366 


1 


530 





2389 





2198 


4 


729(-2) 


4.199 


-5) 


7.341 


-5) 


5.976( 


-5) 





9114 





8600 


-6 


098( 


-2) 





6674 


1 


313 


1 


448 





2515 





2198 


4 


733(-2) 


7.023 


-5) 


1.310 


-4) 


1.182( 


-4) 





8955 





8389 


-6 


954( 


-2) 





6149 


1 


259 


1 


355 





2658 





2199 


4 


745 (-2) 


6.260 


-5) 


1.099 


-4) 


9.640( 


-5) 





8727 





8102 


-8 


066( 


-2) 





5419 


1 


229 


1 


299 





2740 





2199 


4 


756(-2) 


6.800 


-5) 


1.030 


-4) 


9.569( 


-5) 





8560 





7904 


-8 


789 ( 


-2) 





4901 
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TABLE III: Same as Table ini but for diflterent mass binary systems. The prime denotes the values for the more massive star. 



da 


_ 
d 


Q 


M 




Different mass stars, Synchronized case, 7 — 2.5 
J \VE{M)\ \VE{GB)\ \VE{FUS)\ 


02/01 
a'a/a'i 


as/ai 
ag/a'i 




Sec 
Se'c 






X 

x' 












M/R = 


0.14 vs 0.16, 


Mb -- 


= 0.08742 vs. 0.1005 
















2.004 


2.080 


0.1305 


0.1708 


3 


224(-2) 


1.331(-5) 


2.585 


-5) 


2.443( 


-5) 


0.9535 


0.9300 


-1 


939( 


-2) 





8452 
























0.9691 


0.9480 


-1 


740 ( 


-2) 





8833 


1.822 


1.870 


0.1489 


0.1707 


3 


170(-2) 


1.120(-5) 


2.191 


-5) 


2.066( 


-5) 


0.9365 


0.9051 


-2 


617( 


-2) 





7888 
























0.9581 


0.9300 


-2 


318( 


-2) 





8419 


1.640 


1.647 


0.1721 


0.1707 


3 


127(-2) 


8.836(-6) 


1.663 


-5) 


1.558( 


-5) 


0.9072 


0.8651 


-3 


711( 


-2) 





6960 
























0.9398 


0.9017 


-3 


217( 


-2) 





7763 


1.531 


1.500 


0.1894 


0.1706 


3 


113(-2) 


3.528(-5) 


6.227 


-5) 


5.704( 


-5) 


0.8768 


0.8267 


-4 


752 ( 


-2) 





6040 
























0.9218 


0.8758 


-4 


031( 


-2) 





7148 


t 1.482 


1.429 


0.1981 


0.1706 


3 


lll(-2) 


3.086(-5) 


5.309 


-5) 


4.804( 


-5) 


0.8567 


0.8027 


-5 


387( 


-2) 





5443 
























0.9105 


0.8604 


-4 


506( 


-2) 





6775 


1.440 


1.362 


0.2062 


0.1706 


3 


112(-2) 


2.918(-5) 


4.788 


-5) 


4.367( 


-5) 


0.8337 


0.7764 


-6 


056( 


-2) 





4767 
























0.8984 


0.8445 


-4 


989( 


-2) 





6383 


1.421 


1.330 


0.2100 


0.1707 


3 


114(-2) 


3.091(-5) 


4.906 


-5) 


4.524( 


-5) 


0.8209 


0.7624 


-6 


395( 


-2) 





4396 
























0.8921 


0.8365 


-5 


228( 


-2) 





6181 












M/R^ 


0.16 vs 0.18, 


Mb 


= 0.1005 vs. 


1131 
















2.004 


2.179 


0.1364 


0.1917 


3 


904(-2) 


2.247(-5) 


4.390 


-5) 


4.177( 


-5) 


0.9606 


0.9393 


-1 


936( 


-2) 





8634 
























0.9732 


0.9541 


-1 


793 ( 


-2) 





8951 


1.822 


1.965 


0.1554 


0.1916 


3 


839(-2) 


1.908(-5) 


3.759 


-5) 


3.562( 


-5) 


0.9468 


0.9181 


-2 


585( 


-2) 





8150 
























0.9640 


0.9385 


-2 


372( 


-2) 





8588 


1.640 


1.740 


0.1792 


0.1915 


3 


787(-2) 


1.493(-5) 


2.948 


-5) 


2.783( 


-5) 


0.9238 


0.8859 


-3 


608( 


-2) 





7380 
























0.9491 


0.9144 


-3 


256( 


-2) 





8023 


1.495 


1.546 


0.2032 


0.1915 


3 


762(-2) 


2.677(-5) 


4.474 


-5) 


3.698( 


-5) 


0.8908 


0.8413 


-4 


943 ( 


-2) 





6339 
























0.9286 


0.8838 


-4 


357( 


-2) 





7295 


t 1.437 


1.461 


0.2145 


0.1915 


3 


760(-2) 


5.030(-5) 


9.113 


-5) 


8.363( 


-5) 


0.8701 


0.8156 


-5 


715( 


-2) 





5705 
























0.9165 


0.8666 


-4 


967( 


-2) 





6879 


1.386 


1.380 


0.2253 


0.1915 


3 


763(-2) 


4.373(-5) 


7.603 


-5) 


6.901( 


-5) 


0.8447 


0.7858 


-6 


571( 


-2) 





4947 
























0.9025 


0.8479 


-5 


620( 


-2) 





6415 


1.356 


1.329 


0.2330 


0.1915 


3 


768(-2) 


4.533(-5) 


7.403 


-5) 


6.775( 


-5) 


0.8249 


0.7637 


-7 


168( 


-2) 





4364 
























0.8924 


0.8349 


-6 


062( 


-2) 





6086 












M/R^ 


0.18 vs 0.20, 


Mb 


= 0.1131 vs. 


1250 
















1.822 


2.077 


0.1608 


0.2108 


4 


506(-2) 


2.951(-5) 


5.847 


-5) 


5.549( 


-5) 


0.9561 


0.9304 


-2 


583( 


-2) 





8400 
























0.9699 


0.9471 


-2 


465( 


-2) 





8763 


1.640 


1.848 


0.1852 


0.2107 


4 


445(-2) 


4.642(-5) 


8.175 


-5) 


7.062( 


-5) 


0.9381 


0.9029 


-3 


555( 


-2) 





7760 
























0.9578 


0.9268 


-3 


349 ( 


-2) 





8282 


1.459 


1.604 


0.2164 


0.2107 


4 


407(-2) 


3.889(-5) 


6.799 


-5) 


5.682( 


-5) 


0.9050 


0.8567 


-5 


182( 


-2) 





6657 
























0.9366 


0.8935 


-4 


770 ( 


-2) 





7485 


t 1.389 


1.502 


0.2308 


0.2107 


4 


403(-2) 


5.752(-5) 


1.047 


-4) 


9.380( 


-5) 


0.8838 


0.8294 


-6 


127( 


-2) 





5987 
























0.9238 


0.8747 


-5 


560( 


-2) 





7027 


1.349 


1.440 


0.2399 


0.2107 


4 


405(-2) 


5.591(-5) 


9.670 


-5) 


8.621( 


-5) 


0.8672 


0.8091 


-6 


800( 


-21 





5479 
























0.9142 


0.8613 


-6 


114( 


-2) 





6695 


1.313 


1.378 


0.2487 


0.2107 


4 


411(-2) 


5.155(-5) 


8.475 


-5) 


7.485( 


-5) 


0.8476 


0.7860 


-7 


556( 


-2) 





4883 
























0.9035 


0.8468 


-6 


699( 


-2) 





6332 


1.295 


1.345 


0.2534 


0.2107 


4 


415(-2) 


5.715(-5) 


8.631 


-5) 


7.711( 


-5) 


0.8354 


0.7722 


-7 


967( 


-2) 





4522 
























0.8971 


0.8384 


-7 


031( 


-2) 





6120 
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TABLE IV: Same as Table ITTl but for irrotational binaries. 











Identical mass stars, Irrotational case, 7 — 2.5 






da 


d 




M 


J \VE{M)\ \VE{GB)\ \VE{FUS)\ aa/ai aa/ai 




X 



M/R = 0.14 vs 0.14, Mb = 0.08742 vs. 0.08742 



2.185 


2 


266 





1129 





1600 


2 


765(-2) 


1.044 


-5) 


1 


994 ( 


-5) 


1 


865 


-5) 





9678 





9729 


-4 


680 


-4) 





9344 


1.821 


1 


865 





1459 





1598 


2 


599(-2) 


6.162 


-6) 


1 


208( 


-5) 


1 


095 


-5) 





9403 





9474 


-1 


371 


-3) 





8693 


1.565 


1 


552 





1803 





1596 


2 


487(-2) 


1.366 


-5) 


2 


090 ( 


-5) 


1 


565 


-6) 





8903 





9018 


-3 


843 


-3) 





7422 


1.491 


1 


446 





1930 





1595 


2 


458(-2) 


2.819 


-5) 


5 


056( 


-5) 


4 


520 


-5) 





8608 





8751 


-5 


587 


-3) 





6573 


1.454 


1 


383 





1999 





1595 


2 


445 (-2) 


2.273 


-5) 


4 


210( 


-5) 


3 


612 


-5) 





8384 





8548 


-6 


944 


-3) 





5839 


1.436 


1 


346 





2036 





1595 


2 


439(-2) 


1.021 


-5) 


2 


298( 


-5) 


1 


561 


-5) 





8230 





8407 


-7 


880 


-3) 





5286 


1.417 


1 


302 





2075 





1595 


2 


434(-2) 


1.346 


-5) 


1 


706 ( 


-6) 


1 


353 


-5) 





8027 





8219 


-9 


217 


-3) 





4548 


1.406 


1 


272 





2100 





1595 


2 


432(-2) 


4.416 


-5) 


3 


417( 


-5) 


5 


155 


-5) 





7881 





8079 


-1 


036 


-2) 





4062 


















M/R 


= 0.16 


vs 


16, 


Mb 


= 


1005 


vs. 


11005 




















2.003 


2 


159 





1341 





1814 


3 


313(-2) 


1.458 


-5) 


2 


834 ( 


-5) 


2 


643 


-5) 





9630 





9688 


-7 


015 


-4) 





9216 


1.821 


1 


952 





1531 





1813 


3 


215(-2) 


1.132 


-5) 


2 


204( 


-5) 


2 


023 


-5) 





9489 





9557 


-1 


217 


-3) 





8873 


1.638 


1 


735 





1770 





1811 


3 


119(-2) 


6.595 


-6) 


1 


327( 


-5) 


1 


170 


-5) 





9244 





9332 


-2 


361 


-3) 





8262 


1.492 


1 


541 





2014 





1810 


3 


048(-2) 


1.761 


-5) 


2 


864 ( 


-5) 


2 


062 


-5) 





8868 





8992 


-4 


562 


-3) 





7261 


1.455 


1 


486 





2085 





1810 


3 


031(-2) 


4.860 


-5) 


9 


070 ( 


-5) 


8 


259 


-5) 





8719 





8857 


-5 


510 


-3) 





6815 


1.418 


1 


424 





2161 





1809 


3 


016(-2) 


3.861 


-5) 


7 


434 ( 


-5) 


6 


531 


-5) 





8512 





8669 


-6 


816 


-3) 





6118 


1.381 


1 


341 





2242 





1809 


3 


003(-2) 


2.983 


-5) 


1 


013( 


-5) 


3 


139 


-5) 





8160 





8343 


-9 


212 


-3) 





4776 


1.370 


1 


311 





2268 





1809 


2 


999(-2) 


9.006 


-5) 


8 


473 ( 


-5) 


1 


159 


-4) 





8023 





8212 


-1 


048 


-2) 





4336 


















M/R 


= 0.18 


vs 


18, 


Mb 


= 


1131 


vs. 


11131 




















2.003 


2 


268 





1396 





2015 


3 


957(-2) 


2.342 


-5) 


4 


589( 


-5) 


4 


303 


-5) 





9687 





9744 


-6 


386 


-4) 





9336 


1.638 


1 


836 





1837 





2012 


3 


738(-2) 


1.217 


-5) 


2 


452( 


-5) 


2 


199 


-5) 





9377 





9457 


-2 


020 


-3) 





8558 


1.529 


1 


695 





2018 





2011 


3 


676(-2) 


2.773 


-5) 


4 


855( 


-5) 


3 


709 


-5) 





9180 





9278 


-3 


182 


-3) 





8051 


1.456 


1 


595 





2156 





2010 


3 


636(-2) 


2.364 


-5) 


4 


107( 


-5) 


2 


932 


-5) 





8986 





9103 


-4 


452 


-3) 





7522 


1.382 


1 


480 





2313 





2009 


3 


601(-2) 


5.059 


-5) 


9 


859( 


-5) 


8 


570 


-5) 





8674 





8821 


-6 


561 


-3) 





6514 


1.345 


1 


395 





2399 





2009 


3 


583(-2) 


6.626 


-5) 


6 


451( 


-5) 


9 


658 


-5) 





8335 





8505 


-9 


023 


-3) 





5141 


1.327 


1 


347 





2443 





2009 


3 


574(-2) 


2.324 


-4) 


2 


930( 


-4) 


3 


523 


-4) 





8126 





8305 


-1 


139 


-2) 





4510 


1.308 


1 


298 





2489 





2008 


3 


565(-2) 


4.501 


-4) 


5 


946 ( 


-4) 


6 


889 


-4) 





7909 





8094 


-1 


453 


-2) 





4008 


















M/R 


= 0.20 


vs 


20, 


Mb 


= 


1250 


vs. 


D.1250 




















1.821 


2 


178 





1642 





2197 


4 


471(-2) 


3.029 


-5) 


6 


026( 


-5) 


5 


617 


-5) 





9650 





9712 


-9 


649 


-4) 





9219 


1.565 


1 


860 





2011 





2194 


4 


304(-2) 


1.725 


-5) 


3 


517( 


-5) 


3 


136 


-5) 





9407 





9487 


-2 


287 


-3) 





8590 


1.383 


1 


608 





2372 





2192 


4 


196(-2) 


2.926 


-5) 


5 


391( 


-5) 


3 


697 


-5) 





9005 





9124 


-5 


184 


-3) 





7493 


1.346 


1 


549 





2458 





2191 


4 


177(-2) 


2.349 


-5) 


4 


314( 


-5) 


2 


574 


-5) 





8855 





8989 


-6 


292 


-3) 





6993 


1.328 


1 


512 





2504 





2191 


4 


168(-2) 


3.678 


-5) 


8 


046 ( 


-5) 


5 


968 


-5) 





8733 





8877 


-7 


089 


-3) 





6457 


1.309 


1 


464 





2548 





2191 


4 


155(-2) 


9.508 


-5) 


1 


215( 


-4) 


1 


618 


-4) 





8540 





8697 


-8 


684 


-3) 





5621 


1.291 


1 


417 





2593 





2191 


4 


141(-2) 


3.268 


-4) 


4 


764 ( 


-4) 


5 


496 


-4) 





8347 





8515 


-1 


113 


-2) 





5015 


1.272 


1 


370 





2640 





2190 


4 


126(-2) 


6.048 


-4) 


8 


985 ( 


-4) 


1 


013 


-3) 





8157 





8334 


-1 


426 


-2) 





4565 
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TABLE V: Same as Table HVI but for different mass binary systems. The prime denotes the values for the more massive star. 



do 


d 





M 


Different mass stars, Irrotational case, 7 = 2.5 
J |V-E'(M)| \VE{GB)\ \VE{FUS)\ 


0,2 /a\ 
a'2/a'i 


a'z/a'i 


(5ec 
Se'c 




X 

x' 










M/R 


= 0.14 vs 0.16, Mb 


= 0.08742 vs. 


0.1005 












2.003 


2.112 


0.1309 


0.1707 


2.978(-2) 


1.129(-5) 


2.193( 


-5) 


2.036 


-5) 


0.9509 


0.9581 


-1.012( 


-3) 


0.8976 






















0.9677 


0.9725 


-5.307( 


-4) 


0.9307 


1.821 


1.906 


0.1496 


0.1705 


2.888(-2) 


8.442(-6) 


1.665( 


-5) 


1.521 


-5) 


0.9316 


0.9402 


-1.799( 


-3) 


0.8516 






















0.9553 


0.9609 


-9.288( 


-4) 


0.9004 


1.638 


1.687 


0.1732 


0.1704 


2.800(-2) 


4.911(-6) 


9.372( 


-6) 


8.187 


-6) 


0.8967 


0.9083 


-3.597( 


-3) 


0.7634 






















0.9339 


0.9411 


-1.813( 


-3) 


0.8466 


1.528 


1.540 


0.1908 


0.1703 


2.751(-2) 


1.560(-5) 


2.468( 


-5) 


1.816 


-5) 


0.8576 


0.8730 


-6.068( 


-3) 


0.6577 






















0.9116 


0.9209 


-2.942( 


-3) 


0.7884 


1.491 


1.484 


0.1974 


0.1703 


2.735(-2) 


3.031(-5) 


5.431( 


-5) 


4.769 


-5) 


0.8372 


0.8545 


-7.402( 


-3) 


0.5939 






















0.9012 


0.9115 


-3.515( 


-3) 


0.7598 


1.455 


1.420 


0.2045 


0.1702 


2.721(-2) 


1.747(-5) 


3.441( 


-5) 


2.640 


-5) 


0.8071 


0.8270 


-9.393( 


-3) 


0.4892 






















0.8881 


0.8996 


-4.267( 


-3) 


0.7218 


1.436 


1.382 


0.2082 


0.1702 


2.715(-2) 


9.138(-7) 


1.224( 


-5) 


1.942 


-5) 


0.7863 


0.8075 


-1.095( 


-2) 


0.4182 






















0.8798 


0.8922 


-4.740( 


-3) 


0.6966 










M/R 


= 0.16 vs 


18, Mb 


= 


1005 vs. 0.1131 












1.821 


2.002 


0.1562 


0.1913 


3.515(-2) 


1.487(-5) 


2.951( 


-5) 


2.721 


-5) 


0.9428 


0.9508 


-1.527( 


-3) 


0.8751 






















0.9618 


0.9673 


-8.614( 


-4) 


0.9142 


1.638 


1.782 


0.1805 


0.1912 


3.413(-2) 


9.770(-6) 


1.956( 


-5) 


1.748 


-5) 


0.9153 


0.9256 


-2.943( 


-3) 


0.8067 






















0.9442 


0.9510 


-1.622( 


-3) 


0.8697 


1.565 


1.688 


0.1921 


0.1911 


3.374(-2) 


2.498(-5) 


4.281( 


-5) 


3.331 


-5) 


0.8976 


0.9096 


-4.048( 


-3) 


0.7612 






















0.9333 


0.9411 


-2.188( 


-3) 


0.8416 


1.492 


1.588 


0.2051 


0.1910 


3.336(-2) 


2.221(-5) 


3.699( 


-5) 


2.737 


-5) 


0.8726 


0.8870 


-5.704( 


-3) 


0.6917 






















0.9187 


0.9278 


-3.017( 


-3) 


0.8028 


1.455 


1.534 


0.2122 


0.1910 


3.319(-2) 


2.025(-5) 


3.301( 


-5) 


2.328 


-5) 


0.8548 


0.8710 


-6.916( 


-3) 


0.6365 






















0.9092 


0.9192 


-3.593( 


-3) 


0.7766 


1.419 


1.471 


0.2198 


0.1910 


3.302(-2) 


3.680(-5) 


7.472( 


-5) 


6.303 


-5) 


0.8280 


0.8465 


-8.712( 


-3) 


0.5381 






















0.8974 


0.9085 


-4.332( 


-3) 


0.7425 


1.393 


1.417 


0.2254 


0.1909 


3.291(-2) 


2.161(-5) 


6.228( 


-6) 


2.607 


-5) 


0.7996 


0.8197 


-1.111( 


-2) 


0.4382 






















0.8867 


0.8988 


-5.015( 


-3) 


0.7090 










M/R 


= 0.18 vs 


20, Mb 


= 


1131 vs. 0.1250 












1.821 


2.115 


0.1617 


0.2106 


4.146(-2) 


2.407(-5) 


4.781( 


-5) 


4.445 


-5) 


0.9530 


0.9605 


-1.310( 


-3) 


0.8967 






















0.9682 


0.9735 


-7.934( 


-4) 


0.9281 


1.565 


1.798 


0.1984 


0.2103 


3.987(-2) 


3.436(-5) 


6.119( 


-5) 


4.780 


-5) 


0.9182 


0.9287 


-3.278( 


-3) 


0.8090 






















0.9456 


0.9527 


-1.931( 


-3) 


0.8699 


1.492 


1.701 


0.2115 


0.2102 


3.945(-2) 


3.030(-5) 


5.434( 


-5) 


4.047 


-5) 


0.9003 


0.9124 


-4.498( 


-3) 


0.7614 






















0.9345 


0.9426 


-2.605( 


-3) 


0.8405 


1.419 


1.596 


0.2263 


0.2101 


3.904(-2) 


2.503(-5) 


4.469( 


-5) 


3.043 


-5) 


0.8736 


0.8884 


-6.432( 


-3) 


0.6828 






















0.9191 


0.9287 


-3.631( 


-3) 


0.7987 


1.382 


1.534 


0.2343 


0.2101 


3.886(-2) 


1.287(-5) 


2.740( 


-5) 


1.110 


-5) 


0.8507 


0.8675 


-8.001( 


-3) 


0.5961 






















0.9089 


0.9194 


-4.359( 


-3) 


0.7693 


1.364 


1.497 


0.2386 


0.2100 


3.876(-2) 


1.977(-5) 


1.155( 


-5) 


3.520 


-5) 


0.8320 


0.8501 


-9.495( 


-3) 


0.5208 






















0.9026 


0.9138 


-4.810( 


-3) 


0.7503 


1.346 


1.457 


0.2429 


0.2100 


3.866(-2) 


1.021 (-4) 


1.261( 


-4) 


1.631 


-4) 


0.8119 


0.8309 


-1.178( 


-2) 


0.4581 






















0.8951 


0.9070 


-5.342( 


-3) 


0.7258 
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TABLE VI: Same as Table El but for 7 = 2.25. 











Identical mass stars, 


Synchronized case, 7 = 2.25 






da 


d 




M 


J \VE{M)\ 


\VE{GB)\ \VE{FUS)\ az/ai ag/ai 




X 



M/R = 0.12 vs 0.12, Mb = 0.09316 vs. 0.09316 



2 


925 


2 


304 


0.07708 





1734 


3 


724(-2) 


8.099( 


-6) 


1.407( 


-5) 


1.317 


-5) 





9726 





9578 


-1 


402 


-2) 





9002 


2 


701 


2 


113 


0.08635 





1733 


3 


646(-2) 


7.249( 


-6) 


1.259( 


-5) 


1.174 


-5) 





9649 





9457 


-1 


787 


-2) 





8712 


2 


251 


1 


714 


0.1117 





1732 


3 


510(-2) 


4.990( 


-6) 


8.736( 


-6) 


8.072 


-6) 





9354 





9019 


-3 


171 


-2) 





7650 


2 


026 


1 


492 


0.1296 





1731 


3 


462(-2) 


2.062( 


-5) 


3.400( 


-5) 


3.093 


-5) 





9027 





8580 


-4 


536 


-2) 





6549 


1 


936 


1 


392 


0.1383 





1731 


3 


452(-2) 


1.853( 


-5) 


2.870( 


-5) 


2.572 


-5) 





8808 





8307 


-5 


355 


-2) 





5840 


t 1 


881 


1 


327 


0.1440 





1731 


3 


450(-2) 


1.659( 


-5) 


2.505( 


-5) 


2.232 


-5) 





8626 





8091 


-5 


978 


-2) 





5263 


1 


845 


1 


278 


0.1481 





1731 


3 


451(-2) 


1.695( 


-5) 


2.373( 


-5) 


2.153 


-5) 





8471 





7912 


-6 


468 


-2) 





4773 


1 


813 


1 


231 


0.1519 





1731 


3 


455(-2) 


1.961( 


-5) 


2.624( 


-5) 


2.560 


-5) 





8297 





7721 


-6 


961 


-2) 





4231 


















M/R = 


0.14 vs 


0.14, 


Mb -- 


= 0.1080 vs. 


1080 




















2 


701 


2 


229 


0.09129 





1985 


4 


535(-2) 


1.297( 


-5) 


2.327( 


-5) 


2.192 


-5) 





9704 





9535 


-1 


770 


-2) 





8879 


2 


251 


1 


821 


0.1177 





1983 


4 


366(-2) 


9.395( 


-6) 


1.688( 


-5) 


1.573 


-5) 





9468 





9171 


-3 


086 


-2) 





7986 


2 


026 


1 


600 


0.1362 





1982 


4 


302(-2) 


7.380( 


-6) 


1.267( 


-5) 


1.170 


-5) 





9222 





8820 


-4 


332 


-2) 





7109 


1 


936 


1 


505 


0.1451 





1982 


4 


285(-2) 


3.392( 


-5) 


5.874( 


-5) 


5.399 


-5) 





9067 





8611 


-5 


055 


-2) 





6580 


1 


846 


1 


402 


0.1551 





1982 


4 


275(-2) 


2.858( 


-5) 


4.793( 


-5) 


4.339 


-5) 





8849 





8336 


-5 


987 


-2) 





5862 


t 1 


808 


1 


354 


0.1597 





1982 


4 


274(-2) 


2.765( 


-5) 


4.403( 


-5) 


3.958 


-5) 





8726 





8186 


-6 


473 


-2) 





5464 


1 


756 


1 


284 


0.1663 





1982 


4 


277(-2) 


2.930( 


-5) 


4.107( 


-5) 


3.686 


-5) 





8511 





7936 


-7 


250 


-2) 





4779 


1 


724 


1 


235 


0.1707 





1982 


4 


283(-2) 


3.209( 


-5) 


4.785( 


-5) 


4.456 


-5) 





8335 





7740 


-7 


811 


-2) 





4228 


















M/R = 


0.16 vs 


0.16, 


Mb -- 


= 0.1219 vs. 


1219 




















2 


701 


2 


364 


0.09540 





2215 


5 


416(-2) 


2.248( 


-5) 


4.136( 


-5) 


3.929 


-5) 





9754 





9609 


-1 


775 


-2) 





9039 


2 


251 


1 


943 


0.1227 





2213 


5 


214(-2) 


1.703( 


-5) 


3.126( 


-5) 


2.936 


-5) 





9569 





9309 


-3 


044 


-2) 





8295 


2 


026 


1 


720 


0.1417 





2211 


5 


134(-2) 


1.374( 


-5) 


2.462( 


-5) 


2.291 


-5) 





9384 





9029 


-4 


208 


-2) 





7592 


1 


936 


1 


625 


0.1507 





2211 


5 


110(-2) 


2.577( 


-5) 


4.027( 


-5) 


3.303 


-5) 





9271 





8867 


-4 


864 


-2) 





7185 


1 


801 


1 


474 


0.1663 





2211 


5 


086(-2) 


5.052( 


-5) 


9.016( 


-5) 


8.315 


-5) 





9027 





8536 


-6 


189 


-2) 





6334 


t 1 


736 


1 


394 


0.1749 





2211 


5 


082(-2) 


4.393( 


-5) 


7.597( 


-5) 


6.912 


-5) 





8851 





8313 


-7 


040 


-2) 





5746 


1 


666 


1 


298 


0.1850 





2211 


5 


088(-2) 


4.497( 


-5) 


6.625( 


-5) 


5.946 


-5) 





8582 





7992 


-8 


194 


-2) 





4876 


1 


630 


1 


241 


0.1906 





2211 


5 


097(-2) 


5.881( 


-5) 


8.377( 


-5) 


7.738 


-5) 





8384 





7769 


-8 


920 


-2) 





4250 


















M/R = 


0.18 vs 


0.18, 


Mb -- 


= 0.1347 vs. 


1347 




















2 


701 


2 


522 


0.09878 





2419 


6 


251(-2) 


3.691( 


-5) 


6.915( 


-5) 


6.619 


-5) 





9800 





9677 


-1 


825 


-2) 





9190 


2 


251 


2 


084 


0.1267 





2416 


6 


018(-2) 


2.898( 


-5) 


5.451( 


-5) 


5.162 


-5) 





9656 





9434 


-3 


085 


-2) 





8576 


2 


026 


1 


854 


0.1461 





2415 


5 


922(-2) 


2.451( 


-5) 


4.556( 


-5) 


4.254 


-5) 





9517 





9211 


-4 


208 


-2) 





8011 


1 


891 


1 


709 


0.1603 





2414 


5 


878(-2) 


3.901( 


-5) 


6.457( 


-5) 


5.406 


-5) 





9386 





9011 


-5 


190 


-2) 





7505 


1 


712 


1 


502 


0.1832 





2414 


5 


842(-2) 


7.229( 


-5) 


1.321( 


-4) 


1.219 


-4) 





9098 





8605 


-7 


132 


-2) 





6458 


t 1 


662 


1 


439 


0.1904 





2414 


5 


839(-2) 


6.523( 


-5) 


1.173( 


-4) 


1.074 


-4) 





8975 





8445 


-7 


864 


-2) 





6036 


1 


577 


1 


320 


0.2043 





2414 


5 


849(-2) 


5.633( 


-5) 


9.448( 


-5) 


8.565 


-5) 





8673 





8072 


-9 


472 


-2) 





5033 


1 


554 


1 


284 


0.2082 





2414 


5 


855(-2) 


6.519( 


-5) 


1.036( 


-4) 


9.563 


-5) 





8562 





7943 


-9 


980 


-2) 





4677 
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TABLE VII: Same as Table EH but for different mass binary systems. The prime denotes the values for the more massive star. 











Different mass stars, 


Synchronized case, 7 = 2.25 








da 


d 




M 


J \VE{M)\ 


\VE{GB)\ \VE{FUS)\ aa/ai 






X 












aa/ai 


1 1 ' 


5e', 


x' 













MjR = 


0.12 vs 


0.14, 


Mb = 


- 0.09316 vs. 


.1080 




















2.701 


2.169 


0.08889 


0.1859 


4 


063(-2) 


9.821( 


-6) 


1.738( 


-5) 


1.630( 


-5) 





9595 





9400 


-1 


916 


-2) 





8577 





























9744 





9576 


-1 


665 


-2) 





8980 


2.476 


1.970 


0.1005 


0.1858 


3 


983(-2) 


8.511( 


-6) 


1.508( 


-5) 


1.407( 


-5) 





9464 





9208 


-2 


513 


-2) 





8113 





























9663 





9443 


-2 


161 


-2) 





8655 


2.251 


1.763 


0.1148 


0.1858 


3 


912(-2) 


7.046( 


-6) 


1.249( 


-5) 


1.160( 


-5) 





9254 





8914 


-3 


421 


-2) 





7396 





























9538 





9245 


-2 


894 


-2) 





8167 


2.026 


1.539 


0.1330 


0.1857 


3 


858(-2) 


1.413( 


-5) 


1.996( 


-5) 


1.610( 


-5) 





8872 





8419 


-4 


929 


-2) 





6154 





























9322 





8925 


-4 


051 


-2) 





7370 


1.936 


1.440 


0.1418 


0.1857 


3 


845(-2) 


2.219( 


-5) 


3.581( 


-5) 


3.212( 


-5) 





8610 





8105 


-5 


851 


-2) 





5336 





























9187 





8737 


-4 


719 


-2) 





6889 


1.891 


1.385 


0.1466 


0.1857 


3 


842(-2) 


2.220( 


-5) 


3.382( 


-5) 


3.011( 


-5) 





8430 





7899 


-6 


427 


-2) 





4782 





























9100 





8621 


-5 


120 


-2) 





6590 


t 1.866 


1.353 


0.1494 


0.1857 


3 


841(-2) 


2.148( 


-5) 


3.200( 


-5) 


2.882( 


-5) 





8308 





7763 


-6 


788 


-2) 





4410 





























9045 





8549 


-5 


365 


-2) 





6401 












M/R = 


0.14 vs 


0.16, 


Mb -- 


= 0.1080 vs. 


1219 




















2.476 


2.090 


0.1054 


0.2099 


4 


857(-2) 


1.519( 


-5) 


2.769( 


-5) 


2.606( 


-5) 





9564 





9340 


-2 


429 


-2) 





8402 





























9716 





9524 


-2 


183 


-2) 





8829 


2.251 


1.879 


0.1203 


0.2098 


4 


769(-2) 


1.287( 


-5) 


2.347( 


-5) 


2.200( 


-5) 





9404 





9103 


-3 


265 


-2) 





7821 





























9615 





9358 


-2 


899 


-2) 





8415 


2.026 


1.655 


0.1390 


0.2097 


4 


698(-2) 


9.989( 


-6) 


1.804( 


-5) 


1.681( 


-5) 





9128 





8721 


-4 


600 


-2) 





6867 





























9449 





9097 


-4 


003 


-2) 





7761 


1.891 


1.509 


0.1529 


0.2097 


4 


671(-2) 


4.124( 


-5) 


7.177( 


-5) 


6.596( 


-5) 





8840 





8353 


-5 


852 


-2) 





5926 





























9285 





8858 


-4 


988 


-2) 





7155 


1.846 


1.457 


0.1580 


0.2097 


4 


665(-2) 


3.974( 


-5) 


6.562( 


-5) 


5.982( 


-5) 





8704 





8188 


-6 


391 


-2) 





5494 





























9212 





8756 


-5 


397 


-2) 





6893 


t 1.790 


1.386 


0.1649 


0.2097 


4 


662(-2) 


4.089( 


-5) 


5.949( 


-5) 


5.379( 


-5) 





8485 





7932 


-7 


193 


-2) 





4807 





























9101 





8606 


-5 


989 


-2) 





6502 


1.765 


1.353 


0.1682 


0.2097 


4 


662(-2) 


3.651( 


-5) 


5.655( 


-5) 


5.304( 


-5) 





8361 





7793 


-7 


610 


-2) 





4425 





























9043 





8530 


-6 


290 


-2) 





6300 












M/R = 


0.16 vs 


0.18, 


Mb -- 


= 0.1219 vs. 


1347 




















2.251 


2.010 


0.1248 


0.2314 


5 


600(-2) 


2.279( 


-5) 


4.251( 


-5) 


4.018( 


-5) 





9528 





9266 


-3 


174 


-2) 





8189 





























9687 





9465 


-2 


973 


-2) 





8655 


2.026 


1.783 


0.1440 


0.2313 


5 


514(-2) 


3.483( 


-5) 


5.659( 


-5) 


4.793( 


-5) 





9324 





8966 


-4 


399 


-2) 





7440 





























9558 





9253 


-4 


055 


-2) 





8120 


1.891 


1.638 


0.1580 


0.2313 


5 


474(-2) 


3.189( 


-5) 


5.143( 


-5) 


4.267( 


-5) 





9126 





8693 


-5 


496 


-2) 





6747 





























9437 





9065 


-4 


995 


-2) 





7641 


1.801 


1.536 


0.1687 


0.2312 


5 


456(-2) 


3.025( 


-5) 


4.807( 


-5) 


3.943( 


-5) 





8933 





8441 


-6 


480 


-2) 





6098 





























9325 





8899 


-5 


812 


-2) 





7212 


t 1.709 


1.421 


0.1811 


0.2312 


5 


448(-2) 


4.422( 


-5) 


7.306( 


-5) 


6.528( 


-5) 





8638 





8081 


-7 


820 


-2) 





5146 





























9165 





8674 


-6 


886 


-2) 





6626 


1.689 


1.394 


0.1840 


0.2312 


5 


449(-2) 


4.380( 


-5) 


7.051( 


-5) 


6.284( 


-5) 





8551 





7980 


-8 


179 


-2) 





4871 





























9121 





8614 


-7 


165 


-2) 





6468 


1.666 


1.362 


0.1874 


0.2312 


5 


450(-2) 


4.438( 


-5) 


6.860( 


-5) 


6.099( 


-5) 





8440 





7852 


-8 


614 


-2) 





4523 





























9067 





8542 


-7 


496 


-2) 





6277 
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TABLE VIII: Same as Table IVTl but for irrotational binaries. 











Identical mass stars, Irrotational case, 7 = 2.25 






da 


d 




M 


J \VE{M)\ \VE{GB)\ \VE{FUS)\ aa/ai as/ai 




X 



M/R = 0.12 vs 0.12, Mb = 0.09316 vs. 0.09316 



2.924 


2 


331 


0.07722 





1733 


3 


503(-2) 


7 


443(-6) 




1 


278 


-5) 


1.184( 


-5) 





9720 





9764 


-3 


975 


-4) 





9392 


2.699 


2 


145 


0.08654 





1732 


3 


399(-2) 


6 


429(-6) 




1 


096 


-6) 


1.004( 


-5) 





9639 





9689 


-6 


208 


-4) 





9189 


2.249 


1 


756 


0.1121 





1730 


3 


186(-2) 


3 


645(-6) 




5 


956 


-6) 


5.181( 


-6) 





9318 





9390 


-1 


880 


-3) 





8364 


2.023 


1 


540 


0.1302 





1729 


3 


081(-2) 


1 


034(-5) 




1 


390 


-5) 


1.027( 


-5) 





8949 





9051 


-3 


891 


-3) 





7350 


1.933 


1 


440 


0.1390 





1728 


3 


041(-2) 


1 


791(-5) 




2 


913 


-5) 


2.545( 


-5) 





8689 





8813 


-5 


503 


-3) 





6559 


1.887 


1 


384 


0.1438 





1728 


3 


023(-2) 


1 


578(-5) 




2 


575 


-5) 


2.189( 


-5) 





8500 





8640 


-6 


695 


-3) 





5929 


1.842 


1 


315 


0.1490 





1728 


3 


006(-2) 


8 


452(-6) 




1 


842 


-5) 


1.320( 


-5) 





8223 





8385 


-8 


430 


-3) 





4915 


1.818 


1 


272 


0.1518 





1728 


2 


998(-2) 


7 


529(-6) 




4 


168 


-6) 


4.072( 


-6) 





8017 





8194 


-9 


802 


-3) 





4169 
















M/R = 





14 vs 


14, 




Mb 


= 


1080 vs. 


.1080 




















2.699 


2 


260 


0.09153 





1984 


4 


244(-2) 


1 


166(-5) 




2 


069 


-5) 


1.923( 


-5) 





9696 





9746 


-5 


469 


-4) 





9320 


2.249 


1 


863 


0.1182 





1981 


3 


989(-2) 


7 


237(-6) 




1 


255 


-5) 


1.123( 


-5) 





9440 





9507 


-1 


557 


-3) 





8653 


1.933 


1 


555 


0.1460 





1979 


3 


813(-2) 


1 


479(-5) 




2 


111 


-5) 


1.534( 


-5) 





8979 





9084 


-4 


220 


-3) 





7377 


1.843 


1 


451 


0.1562 





1978 


3 


767(-2) 


2 


958(-5) 




5 


118 


-5) 


4.526( 


-5) 





8710 





8838 


-6 


024 


-3) 





6533 


1.797 


1 


390 


0.1617 





1978 


3 


745(-2) 


2 


347(-5) 




4 


246 


-5) 


3.576( 


-5) 





8502 





8647 


-7 


406 


-3) 





5789 


1.774 


1 


353 


0.1647 





1977 


3 


735(-2) 


1 


372(-5) 




3 


164 


-5) 


2.314( 


-5) 





8352 





8509 


-8 


379 


-3) 





5209 


1.752 


1 


309 


0.1677 





1977 


3 


726(-2) 


1 


877(-5) 




2 


257 


-6) 


1.639( 


-5) 





8148 





8320 


-9 


855 


-3) 





4447 


1.742 


1 


289 


0.1690 





1977 


3 


722(-2) 


4 


311(-5) 




2 


740 


-5) 


4.625( 


-5) 





8052 





8229 


-1 


068 


-2) 





4127 
















M/R = 





16 vs 


16, 




Mb 


= 


1219 vs. 


.1219 




















2.699 


2 


396 


0.09567 





2213 


5 


089(-2) 


1 


986(-5) 




3 


624 


-5) 


3.404( 


-5) 





9749 





9799 


-4 


866 


-4) 





9442 


2.249 


1 


985 


0.1232 





2210 


4 


795(-2) 


1 


381(-5) 




2 


503 


-5) 


2.292( 


-5) 





9548 





9611 


-1 


302 


-3) 





8911 


2.024 


1 


769 


0.1425 





2208 


4 


649(-2) 


9 


182(-6) 




1 


624 


-5) 


1.431( 


-5) 





9339 





9418 


-2 


449 


-3) 





8342 


1.934 


1 


677 


0.1518 





2207 


4 


592(-2) 


2 


305(-5) 




3 


618 


-5) 


2.759( 


-5) 





9208 





9297 


-3 


308 


-3) 





7978 


1.843 


1 


580 


0.1620 





2206 


4 


536(-2) 


2 


066(-5) 




3 


166 


-5) 


2.288( 


-5) 





9030 





9134 


-4 


559 


-3) 





7458 


1.753 


1 


472 


0.1737 





2206 


4 


483(-2) 


4 


918(-5) 




8 


947 


-5) 


8.036( 


-5) 





8757 





8885 


-6 


554 


-3) 





6574 


1.707 


1 


404 


0.1800 





2205 


4 


459(-2) 


2 


734(-5) 




5 


954 


-5) 


4.726( 


-5) 





8515 





8664 


-8 


228 


-3) 





5623 


1.662 


1 


307 


0.1868 





2205 


4 


434(-2) 


1 


526(-4) 




1 


638 


-4) 


2.070( 


-4) 





8080 





8257 


-1 


237 


-2) 





4113 
















M/R = 





18 vs 


18, 




Mb 


= 


1347 vs. 


.1347 




















2.699 


2 


554 


0.09906 





2417 


5 


897(-2) 


3 


207(-5) 




5 


974 


-5) 


5.656( 


-5) 





9798 





9846 


-4 


427 


-4) 





9555 


2.249 


2 


126 


0.1274 





2413 


5 


570(-2) 


2 


405(-5) 




4 


480 


-5) 


4.156( 


-5) 





9643 





9702 


-1 


115 


-3) 





9137 


1.889 


1 


765 


0.1614 





2410 


5 


310(-2) 


1 


304(-5) 




2 


334 


-5) 


2.046( 


-5) 





9337 





9419 


-3 


042 


-3) 





8284 


1.708 


1 


562 


0.1849 





2408 


5 


189(-2) 


2 


551(-5) 




4 


151 


-5) 


2.837( 


-5) 





8982 





9095 


-5 


944 


-3) 





7225 


1.663 


1 


502 


0.1918 





2407 


5 


161(-2) 


5 


802(-5) 




1 


087 


-4) 


9.491( 


-5) 





8826 





8952 


-7 


211 


-3) 





6665 


1.640 


1 


467 


0.1953 





2407 


5 


147(-2) 


4 


481(-5) 




8 


930 


-5) 


7.347( 


-5) 





8704 





8841 


-8 


092 


-3) 





6150 


1.618 


1 


420 


0.1989 





2407 


5 


132(-2) 


4 


391(-5) 




3 


902 


-5) 


6.787( 


-5) 





8509 





8661 


-9 


735 


-3) 





5322 


1.595 


1 


369 


0.2026 





2406 


5 


113(-2) 


2 


193(-4) 




2 


871 


-4) 


3.426( 


-4) 





8288 





8456 


-1 


248 


-2) 





4615 
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TABLE IX: Same as Table IVTTTl but for different mass binary systems. Tfie prime denotes tfie values for the more massive star. 



da 


d 


n 


M 


Different mass stars, Irrotational case, 7 — 2.25 
J \VE{M)\ \VE{GB)\ \VE{FUS)\ 


a-i/ax 
aa/a'i 


as/ai 
ag/a'i 


5ec 
&e!c 




X 

x' 










M/R -- 


= 0.12 vs 0.14, Mb -- 


= 0.09316 vs. 0.1080 












2.699 


2.200 


0.08910 


0.1858 


3.794(-2) 


8.995(-6) 1.576 


-5) 


1.458( 


-5) 


0.9582 


0.9643 


-8.350( 


-4) 


0.9072 




















0.9737 


0.9778 


-4.042( 


-4) 


0.9404 


2.474 


2.006 


0.1008 


0.1857 


3.678(-2) 


7.210(-6) 1.246 


-5) 


1.134( 


-5) 


0.9440 


0.9512 


-1.386( 


-3) 


0.8716 




















0.9651 


0.9697 


-6.597( 


-4) 


0.9182 


2.249 


1.805 


0.1152 


0.1856 


3.561(-2) 


5.310(-6) 8.974 


-6) 


7.943( 


-6) 


0.9210 


0.9298 


-2.511( 


-3) 


0.8119 




















0.9515 


0.9570 


-1.166( 


-3) 


0.8822 


2.023 


1.587 


0.1337 


0.1854 


3.447(-2) 


1.322(-5) 1.850 


-5) 


1.394( 


-5) 


0.8772 


0.8896 


-5.215( 


-3) 


0.6903 




















0.9273 


0.9346 


-2.322( 


-3) 


0.8165 


1.933 


1.487 


0.1427 


0.1854 


3.403(-2) 


2.179(-5) 3.612 


-5) 


3.144( 


-5) 


0.8447 


0.8598 


-7.435( 


-3) 


0.5874 




















0.9117 


0.9203 


-3.193( 


-3) 


0.7718 


1.887 


1.429 


0.1476 


0.1853 


3.382(-2) 


1.794(-5) 3.037 


-5) 


2.523( 


-5) 


0.8190 


0.8363 


-9.171( 


-3) 


0.4967 




















0.9014 


0.9108 


-3.804( 


-3) 


0.7411 


1.865 


1.395 


0.1501 


0.1853 


3.373(-2) 


1.079(-5) 2.103 


-5) 


1.498( 


-5) 


0.8010 


0.8196 


-1.042( 


-2) 


0.4320 




















0.8953 


0.9052 


-4.177( 


-3) 


0.7220 










M/R 


= 0.14 vs 0.16, Mb 


= 


1080 vs. 


.1219 












2.474 


2.126 


0.1058 


0.2097 


4.508(-2) 


1.322(-5) 2.378 


-5) 


2.199( 


-5) 


0.9546 


0.9613 


-1.128( 


-3) 


0.8959 




















0.9707 


0.9753 


-5.989( 


-4) 


0.9314 


2.249 


1.921 


0.1208 


0.2096 


4.371(-2) 


1.029(-5) 1.834 


-5) 


1.665( 


-5) 


0.9370 


0.9449 


-1.971( 


-3) 


0.8498 




















0.9598 


0.9651 


-1.026( 


-3) 


0.9020 


2.024 


1.704 


0.1398 


0.2094 


4.235(-2) 


6.653(-6) 1.100 


-5) 


9.575( 


-6) 


0.9055 


0.9158 


-3.843( 


-3) 


0.7620 




















0.9411 


0.9477 


-1.950( 


-3) 


0.8509 


1.888 


1.560 


0.1539 


0.2093 


4.156(-2) 


1.779(-5) 2.613 


-5) 


1.900( 


-5) 


0.8706 


0.8841 


-6.333( 


-3) 


0.6606 




















0.9221 


0.9303 


-3.088( 


-3) 


0.7974 


1.843 


1.506 


0.1592 


0.2092 


4.132(-2) 


3.381(-5) 5.868 


-5) 


5.142( 


-5) 


0.8523 


0.8674 


-7.640( 


-3) 


0.5990 




















0.9135 


0.9224 


-3.646( 


-3) 


0.7720 


1.797 


1.444 


0.1648 


0.2092 


4.108(-2) 


2.288(-5) 4.357 


-5) 


3.482( 


-5) 


0.8245 


0.8418 


-9.591( 


-3) 


0.4945 




















0.9028 


0.9126 


-4.365( 


-3) 


0.7394 


1.775 


1.407 


0.1678 


0.2092 


4.097(-2) 


1.923(-6) 1.354 


-5) 


1.255( 


-6) 


0.8037 


0.8225 


-1.127( 


-2) 


0.4210 




















0.8963 


0.9066 


-4.813( 


-3) 


0.7184 










M/R 


= 0.16 vs 0.18, Mb 


= 


1219 vs. 


.1347 












2.249 


2.053 


0.1253 


0.2312 


5.166(-2) 


1.895(-5) 3.502 


-5) 


3.236( 


-5) 


0.9503 


0.9576 


-1.579( 


-3) 


0.8813 




















0.9675 


0.9727 


-9.161( 


-4) 


0.9207 


2.024 


1.833 


0.1448 


0.2310 


5.011(-2) 


1.339(-5) 2.445 


-5) 


2.196( 


-5) 


0.9274 


0.9364 


-2.956( 


-3) 


0.8192 




















0.9533 


0.9594 


-1.667( 


-3) 


0.8812 


1.934 


1.740 


0.1542 


0.2309 


4.951(-2) 


2.925(-5) 4.763 


-5) 


3.735( 


-5) 


0.9129 


0.9231 


-3.985( 


-3) 


0.7790 




















0.9445 


0.9514 


-2.218( 


-3) 


0.8570 


1.798 


1.591 


0.1702 


0.2307 


4.863(-2) 


2.434(-5) 3.878 


-5) 


2.804( 


-5) 


0.8797 


0.8927 


-6.541( 


-3) 


0.6784 




















0.9261 


0.9344 


-3.513( 


-3) 


0.8038 


1.753 


1.535 


0.1762 


0.2307 


4.835(-2) 


2.135(-5) 3.350 


-5) 


2.248( 


-5) 


0.8614 


0.8761 


-7.923( 


-3) 


0.6138 




















0.9175 


0.9265 


-4.165( 


-3) 


0.7783 


1.708 


1.467 


0.1826 


0.2306 


4.808(-2) 


8.588(-7) 1.681 


-5) 


4.205( 


-7) 


0.8290 


0.8463 


-1.042( 


-2) 


0.4858 




















0.9068 


0.9167 


-5.008( 


-3) 


0.7446 


1.694 


1.443 


0.1845 


0.2306 


4.800(-2) 


3.132(-5) 2.533 


-5) 


4.673( 


-5) 


0.8158 


0.8340 


-1.178( 


-2) 


0.4422 




















0.9028 


0.9131 


-5.317( 


-3) 


0.7318 
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TABLE X: Same as TableEllbut for 7 = 1.8. 

Identical mass stars, Synchronized case, 7 = 1.8 



da 




d 




ft 




M 




J 


\VE{M)\ 


\VE{GB)\ 


\VE{FUS)\ 


ai/a\ 


aa/ai 




(5ec 






X 


















M/R 


= 0.08 


vs 


08, 


Mb 


= 


1485 vs. 


0.1485 




















7.747 


2 


340 





02361 





2852 





1200 


3.257( 


-6) 


4 


951( 


-6) 


4.521 


;-6) 





9771 





9657 


-1 


777 


-2) 





9005 


6.887 


2 


058 





02800 





2850 





1155 


2.891( 


-6) 


4 


363( 


-6) 


3.941 


'-6) 





9668 





9499 


-2 


546 


-2) 





8546 


6.026 


1 


763 





03396 





2849 





1112 


2.427( 


-6) 


3 


637( 


-6) 


3.258 


;-6) 





9475 





9218 


-3 


873 


-2) 





7729 


5.596 


1 


605 





03779 





2848 





1092 


2.044( 


-6) 


3 


034( 


-6) 


2.689 


;-6) 





9305 





8985 


-4 


931 


-2) 





7047 


5.165 


1 


434 





04243 





2848 





1074 


3.124( 


-6) 


5 


190( 


-6) 


4.906 


'-6) 





9026 





8625 


-6 


482 


-2) 





5977 


4.907 


1 


316 





04571 





2847 





1066 


2.566( 


-6) 


3 


873( 


-6) 


3.543 


;-6) 





8743 





8285 


-7 


806 


-2) 





4951 


4.820 


1 


271 





04692 





2847 





1064 


7.892( 


-6) 


9 


143( 


-6) 


7.958 


;-6) 





8604 





8128 


-8 


351 


-2) 





4476 


4.734 


1 


222 





04818 





2847 





1062 


1.199( 


-5) 


1 


547( 


-5) 


1.478 


;-5) 





8432 





7937 


-8 


954 


-2) 





3890 


















M/R 


= 0.10 


vs 


10, 


Mb 


= 


1725 vs. 


0.1725 




















6.887 


2 


319 





02980 





3281 





1432 


6.430( 


-6) 


1 


044( 


-5) 


9.675 


;-6) 





9768 





9649 


-2 


239 


-2) 





8962 


6.026 


2 


004 





03608 





3279 





1376 


5.553( 


-6) 


8 


969( 


-6) 


8.212 


;-6) 





9646 





9459 


-3 


351 


-2) 





8404 


5.596 


1 


840 





04010 





3278 





1349 


4.987( 


-6) 


8 


040 ( 


-6) 


7.310 


;-6) 





9544 





9309 


-4 


212 


-2) 





7959 


5.166 


1 


668 





04494 





3277 





1323 


4.330( 


-6) 


7 


046 ( 


-6) 


6.344 


'-6) 





9390 





9089 


-5 


424 


-2) 





7311 


4.908 


1 


559 





04834 





3277 





1310 


9.463( 


-6) 


1 


178( 


-5) 


9.255 


;-6) 





9253 





8903 


-6 


412 


-2) 





6760 


4.563 


1 


401 





05362 





3276 





1294 


1.203( 


-5) 


1 


661( 


-5) 


1.410 


;-5) 





8970 





8541 


-8 


209 


-2) 





5678 


4.391 


1 


311 





05665 





3275 





1288 


1.290( 


-5) 


1 


597( 


-5) 


1.336 


'-5) 





8743 





8269 


-9 


420 


-2) 





4855 


4.305 


1 


260 





05828 





3275 





1286 


1.476( 


-5) 


2 


046 ( 


-5) 


1.787 


;-5) 





8586 





8090 


-1 


013 


-1) 





4309 


















M/R 


= 0.12 


vs 


12, 


Mb 


= 


1920 vs. 


0.1920 




















6.887 


2 


607 





03112 





3622 





1666 


1.307( 


-5) 


2 


246 ( 


-5) 


2.123 


;-5) 





9835 





9753 


-2 


125 


-2) 





9255 


6.026 


2 


264 





03764 





3620 





1598 


1.153( 


-5) 


1 


970( 


-5) 


1.841 


;-5) 





9755 





9623 


-3 


144 


-2) 





8865 


5.166 


1 


907 





04679 





3617 





1533 


9.487( 


-6) 


1 


618( 


-5) 


1.490 


'-5) 





9597 





9377 


-4 


980 


-2) 





8129 


4.736 


1 


718 





05287 





3616 





1503 


8.299( 


-6) 


1 


421( 


-5) 


1.296 


'-5) 





9451 





9162 


-6 


503 


-2) 





7487 


4.306 


1 


514 





06043 





3614 





1478 


1.520( 


-5) 


2 


794( 


-5) 


2.682 


;-5) 





9198 





8813 


-8 


808 


-2) 





6444 


4.047 


1 


375 





06589 





3614 





1466 


2.016( 


-5) 


3 


009( 


-5) 


2.594 


'-5) 





8928 





8471 


-1 


084 


-1) 





5414 


3.961 


1 


324 





06791 





3613 





1464 


2.084( 


-5) 


3 


036( 


-5) 


2.618 


;-5) 





8797 





8314 


-1 


169 


-1) 





4939 


3.875 


1 


268 





07004 





3613 





1461 


2.752( 


-5) 


4 


055( 


-5) 


3.707 


;-5) 





8630 





8120 


-1 


262 


-1) 





4351 


















M/R 


= 0.14 


vs 


14, 


Mb 


= 


2065 vs. 


0.2065 




















6.456 


2 


748 





03511 





3868 





1803 


2.378( 


-5) 


4 


247( 


-5) 


4.065 


'-5) 





9860 





9791 


-2 


779 


-2) 





9353 


5.165 


2 


168 





04808 





3863 





1690 


1.918( 


-5) 


3 


423( 


-5) 


3.209 


;-5) 





9728 





9572 


-5 


268 


-2) 





8686 


4.305 


1 


754 





06196 





3860 





1622 


2.654( 


-5) 


3 


947( 


-5) 


3.317 


'-5) 





9493 





9209 


-8 


992 


-2) 





7588 


3.875 


1 


523 





07163 





3858 





1595 


2.584( 


-5) 


3 


656( 


-5) 


2.994 


;-5) 





9225 





8832 


-1 


238 


-1) 





6452 


3.703 


1 


419 





07625 





3857 





1587 


3.486( 


-5) 


5 


796 ( 


-5) 


5.063 


;-5) 





9042 





8592 


-1 


428 


-1) 





5729 


3.617 


1 


363 





07876 





3857 





1584 


3.603( 


-5) 


5 


962( 


-5) 


5.219 


;-5) 





8915 





8436 


-1 


541 


-1) 





5259 


3.531 


1 


302 





08141 





3857 





1582 


4.031( 


-5) 


6 


668( 


-5) 


5.910 


;-5) 





8753 





8243 


-1 


668 


-1) 





4674 


3.488 


1 


269 





08280 





3857 





1581 


5.067( 


-5) 


8 


384( 


-5) 


7.517 


;-5) 





8653 





8127 


-1 


735 


-1) 





4325 
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TABLE XI: Same as Table 1X1 but for different mass binary systems. The prime denotes the values for the more massive star. 



do 


d 




M 




Different mass stars, Synchronized case, 7 = 1.8 
J \yE{M)\ \VE{GB)\ \VE{FUS)\ 


02/01 

aa/a'i 


aa/ai 
a's/a'i 




Sec 
Se'c 






X 

x' 












M/R 


= 0.08 


vs 


10, 


Mb 


= 


1485 vs. 


0.1725 
















7.747 


2.473 


0.02440 


0.3067 





1335 


5.001 


-6) 


7 


956( 


-6) 


7.382 


'-6) 


0.9734 


0.9617 


-1 


912 


-2) 





8890 




























0.9860 


0.9781 


-1 


478 


-2) 





9354 


6.887 


2.178 


0.02892 


0.3066 





1285 


4.532 


-6) 


7 


170( 


-6) 


6.587 


;-6) 


0.9614 


0.9440 


-2 


745 


-2) 





8377 




























0.9801 


0.9683 


-2 


097 


-2) 





9064 


6.457 


2.027 


0.03174 


0.3065 





1259 


4.226 


-6) 


6 


671( 


-6) 


6.095 


;-6) 


0.9521 


0.9307 


-3 


360 


-2) 





7990 




























0.9756 


0.9610 


-2 


545 


-2) 





8849 


6.026 


1.872 


0.03505 


0.3064 





1235 


3.877 


-6) 


6 


127( 


-6) 


5.562 


;-6) 


0.9388 


0.9124 


-4 


190 


-2) 





7457 




























0.9695 


0.9513 


-3 


136 


-2) 





8561 


5.596 


1.709 


0.03898 


0.3063 





1212 


3.464 


-6) 


5 


487( 


-6) 


4.950 


;-6) 


0.9188 


0.8858 


-5 


353 


-2) 





6682 




























0.9607 


0.9377 


-3 


937 


-2) 





8161 


5.166 


1.533 


0.04373 


0.3063 





1191 


3.302 


-6) 


5 


186( 


-6) 


4.666 


;-6) 


0.8850 


0.8439 


-7 


073 


-2) 





5442 




























0.9474 


0.9180 


-5 


064 


-2) 





7576 


4.907 


1.412 


0.04708 


0.3062 





1181 


1.115 


-5) 


1 


408 ( 


-5) 


1.234 


;-5) 


0.8485 


0.8021 


-8 


565 


-2) 





4191 




























0.9355 


0.9012 


-5 


985 


-2) 





7077 












M/R 


= 0.10 


vs 


12, 


Mb 


= 


1725 vs. 


0.1920 
















6.887 


2.454 


0.03047 


0.3452 





1543 


9.525 


-6) 


1 


601( 


-5) 


1.504 


;-5) 


0.9742 


0.9621 


-2 


353 


-2) 





8881 




























0.9852 


0.9770 


-2 


031 


-2) 





9306 


6.026 


2.124 


0.03687 


0.3450 





1482 


8.344 


-6) 


1 


394( 


-5) 


1.295 


;-5) 


0.9606 


0.9417 


-3 


526 


-2) 





8279 




























0.9780 


0.9649 


-3 


005 


-2) 





8943 


5.166 


1.777 


0.04589 


0.3447 





1423 


6.858 


-6) 


1 


155( 


-5) 


1.056 


;-5) 


0.9321 


0.9015 


-5 


722 


-2) 





7095 




























0.9637 


0.9420 


-4 


757 


-2) 





8257 


4.736 


1.587 


0.05189 


0.3446 





1398 


6.736 


-6) 


1 


128( 


-5) 


1.016 


;-5) 


0.9030 


0.8636 


-7 


636 


-2) 





5973 




























0.9505 


0.9220 


-6 


204 


-2) 





7658 


4.563 


1.504 


0.05469 


0.3446 





1389 


1.792 


-5) 


2 


685( 


-5) 


2.361 


'-5) 


0.8844 


0.8408 


-8 


701 


-2) 





5296 




























0.9429 


0.9110 


-6 


981 


-2) 





7329 


4.477 


1.460 


0.05618 


0.3445 





1385 


1.819 


-5) 


2 


616( 


-5) 


2.282 


'-5) 


0.8724 


0.8267 


-9 


316 


-2) 





4873 




























0.9384 


0.9046 


-7 


415 


-2) 





7137 


4.391 


1.414 


0.05775 


0.3445 





1381 


2.022 


-5) 


2 


732 ( 


-5) 


2.379 


;-5) 


0.8577 


0.8099 


-9 


998 


-2) 





4366 




























0.9332 


0.8974 


-7 


888 


-2) 





6922 












M/R 


= 0.12 


vs 


14, 


Mb 


= 


1920 vs. 


0.2065 
















5.596 


2.217 


0.04239 


0.3742 





1643 


1.547 


-5) 


2 


693( 


-5) 


2.543 


;-5) 


0.9668 


0.9498 


-4 


044 


-2) 





8488 




























0.9803 


0.9685 


-4 


070 


-2) 





9029 


5.166 


2.028 


0.04744 


0.3740 





1609 


1.408 


-5) 


2 


449 ( 


-5) 


2.298 


;-5) 


0.9568 


0.9346 


-5 


145 


-2) 





8035 




























0.9747 


0.9592 


-5 


124 


-2) 





8746 


4.736 


1.832 


0.05359 


0.3739 





1576 


1.263 


-5) 


2 


166( 


-5) 


2.077 


'-5) 


0.9411 


0.9119 


-6 


719 


-2) 





7359 




























0.9663 


0.9456 


-6 


600 


-2) 





8333 


4.305 


1.622 


0.06120 


0.3737 





1548 


2.080 


-5) 


2 


979( 


-5) 


2.463 


;-5) 


0.9137 


0.8748 


-9 


125 


-2) 





6256 




























0.9527 


0.9245 


-8 


739 


-2) 





7696 


4.047 


1.483 


0.06670 


0.3736 





1534 


2.426 


-5) 


3 


658( 


-5) 


3.140 


;-5) 


0.8843 


0.8383 


-1 


124 


-1) 





5161 




























0.9397 


0.9055 


-1 


053 


-1) 





7123 


3.961 


1.432 


0.06872 


0.3736 





1530 


2.527 


-5) 


3 


762 ( 


-5) 


3.233 


'-5) 


0.8700 


0.8213 


-1 


211 


-1) 





4652 




























0.9340 


0.8975 


-1 


125 


-1) 





6880 


3.875 


1.378 


0.07086 


0.3736 





1526 


3.746 


-5) 


4 


675( 


-5) 


4.142 


;-5) 


0.8518 


0.8007 


-1 


300 


-1) 





4028 




























0.9273 


0.8883 


-1 


205 


-1) 





6601 
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TABLE XII: Same as Table|X|but for irrotational binaries. 











Identical mass stars, Irrotational case, 7 = 1.8 






da 


d 




M 


J \VE{M)\ \VE{GB)\ \VE{FUS)\ aa/ai ag/ai 


5ec 


X 



















M/R 


= 0.08 


vs 


08, 


Mb 


= 


1485 


vs. 


0.1485 




















7.746 


2 


369 





02363 





2851 


0.1147 


3.139 


-6) 


4 


714( 


-6) 


4 


269 


;-6) 





9772 





9804 


-3 


530 


-4) 





9387 


6.024 


1 


812 





03401 





2848 


0.1033 


2.178 


-6) 


3 


112( 


-6) 


2 


689 


;-6) 





9480 





9527 


-1 


476 


-3) 





8488 


5.593 


1 


663 





03785 





2847 


0.1004 


1.839 


-6) 


2 


535( 


-6) 


2 


128 


;-6) 





9317 





9373 


-2 


346 


-3) 





7966 


5.162 


1 


504 





04251 





2845 


0.09736 


7.196 


-6) 


9 


673 ( 


-6) 


8 


089 


;-6) 





9050 





9122 


-4 


010 


-3) 





7077 


4.989 


1 


433 





04465 





2845 


0.09616 


7.030 


-6) 


9 


188( 


-6) 


7 


585 


;-6) 





8887 





8968 


-5 


091 


-3) 





6494 


4.816 


1 


354 





04699 





2844 


0.09497 


7.099 


-6) 


9 


516( 


-6) 


7 


831 


;-6) 





8648 





8745 


-6 


631 


-3) 





5585 


4.730 


1 


306 





04825 





2844 


0.09439 


6.175 


-6) 


9 


979 ( 


-6) 


8 


045 


;-6) 





8468 





8576 


-7 


741 


-3) 





4858 


4.643 


1 


245 





04959 





2844 


0.09384 


5.887 


-7) 


1 


112( 


-5) 


7 


910 


;-6) 





8190 





8315 


-9 


440 


-3) 





3770 


















M/R 


= 0.10 


vs 


10, 


A?B 


= 


1725 


vs. 


0.1725 




















6.885 


2 


353 





02983 





3280 





1365 


6.157 


-6) 


9 


913( 


-6) 


9 


110 


;-6) 





9769 





9806 


-4 


665 


-4) 





9372 


5.163 


1 


730 





04504 





3275 





1221 


3.685 


-6) 


5 


585( 


-6) 


4 


817 


;-6) 





9401 





9457 


-2 


309 


-3) 





8202 


4.732 


1 


557 





05100 





3273 





1183 


3.314 


-6) 


4 


566( 


-6) 


3 


868 


;-6) 





9155 





9225 


-3 


970 


-3) 





7353 


4.559 


1 


483 





05377 





3272 





1168 


1.302 


-5) 


1 


906 ( 


-5) 


1 


611 


;-5) 





9004 





9085 


-5 


156 


-3) 





6850 


4.387 


1 


400 





05682 





3271 





1153 


1.205 


-5) 


1 


747( 


-5) 


1 


439 


;-5) 





8787 





8881 


-6 


775 


-3) 





6030 


4.300 


1 


352 





05846 





3271 





1146 


1.163 


-5) 


1 


773 ( 


-5) 


1 


438 


;-5) 





8626 





8730 


-7 


880 


-3) 





5375 


4.257 


1 


324 





05931 





3271 





1142 


1.048 


-5) 


1 


804 ( 


-5) 


1 


426 


;-5) 





8516 





8628 


-8 


589 


-3) 





4921 


4.213 


1 


291 





06019 





3271 





1139 


3.464 


-6) 


1 


495 ( 


-5) 


9 


870 


;-6) 





8371 





8492 


-9 


611 


-3) 





4333 


















M/R 


= 0.12 


vs 


12, 


Mb 


= 


1920 


vs. 


0.1920 




















6.885 


2 


640 





03115 





3621 





1600 


1.266 


-5) 


2 


172( 


-5) 


2 


045 


;-5) 





9837 





9874 


-3 


748 


-4) 





9574 


5.163 


1 


965 





04691 





3615 





1434 


8.636 


-6) 


1 


438( 


-5) 


1 


300 


;-5) 





9604 





9654 


-1 


553 


-3) 





8831 


4.733 


1 


787 





05303 





3613 





1390 


7.042 


-6) 


1 


150( 


-5) 


1 


015 


;-5) 





9464 





9522 


-2 


537 


-3) 





8373 


4.302 


1 


599 





06064 





3610 





1346 


5.564 


-6) 


8 


300 ( 


-6) 


7 


114 


;-6) 





9227 





9297 


-4 


481 


-3) 





7542 


4.043 


1 


475 





06617 





3608 





1320 


2.167 


-5) 


3 


418( 


-5) 


2 


909 


;-5) 





8983 





9069 


-6 


804 


-3) 





6687 


3.870 


1 


375 





07034 





3607 





1303 


1.597 


-5) 


2 


778 ( 


-5) 


2 


177 


;-5) 





8685 





8790 


-9 


350 


-3) 





5475 


3.827 


1 


342 





07146 





3607 





1298 


1.084 


-5) 


2 


310( 


-5) 


1 


594 


;-5) 





8555 





8669 


-1 


028 


-2) 





4923 


3.783 


1 


302 





07260 





3607 





1294 


1.727 


-5) 


4 


221( 


-7) 


1 


228 


;-5) 





8366 





8495 


-1 


204 


-2) 





4192 


















M/R 


= 0.14 


vs 


14, 


Mb 


= 


2065 


vs. 


0.2065 




















6.025 


2 


602 





03877 





3865 





1691 


2.183 


-5) 


3 


889( 


-5) 


3 


688 


;-5) 





9833 





9874 


-6 


324 


-4) 





9554 


4.302 


1 


840 





06217 





3856 





1501 


1.315 


-5) 


2 


284( 


-5) 


2 


060 


;-5) 





9518 





9576 


-3 


291 


-3) 





8510 


3.872 


1 


633 





07197 





3853 





1452 


2.359 


-5) 


3 


504 ( 


-5) 


2 


722 


;-5) 





9282 





9354 


-6 


118 


-3) 





7701 


3.699 


1 


543 





07664 





3851 





1433 


2.193 


-5) 


3 


218( 


-5) 


2 


409 


;-5) 





9128 





9210 


-8 


110 


-3) 





7146 


3.613 


1 


495 





07919 





3850 





1423 


3.590 


-5) 


6 


044 ( 


-5) 


5 


216 


;-5) 





9025 





9113 


-9 


426 


-3) 





6751 


3.527 


1 


441 





08186 





3850 





1414 


3.051 


-5) 


5 


176( 


-5) 


4 


295 


;-5) 





8884 





8981 


-1 


111 


-2) 





6173 


3.440 


1 


372 





08469 





3849 





1404 


2.422 


-6) 


1 


779 ( 


-5) 


4 


755 


;-5) 





8629 





8744 


-1 


382 


-2) 





5056 


3.414 


1 


344 





08555 





3849 





1401 


3.859 


-5) 


2 


961( 


-5) 


4 


862 


;-5) 





8498 





8624 


-1 


561 


-2) 





4548 
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TABLE XIII: Same as Table PTm but for different mass binary systems. The prime denotes the values of the more massive star. 



da 


d 




M 




Different mass stars, Irrotational case, 7 — 1.8 
J \VE{M)\ \VE{GB)\ \VE{FUS)\ 


ai/ou\ 
a'2/a.'i 


as/ 0,1 
a'z/a'i 




5ec 
5e'c 






X 

x' 












M/R 


= 0.08 


vs 


10, 


Mb 


= 


1485 vs. 


0.1725 
















6.885 


2.215 


0.02895 


0.3065 





1219 


4.305 


-6) 


6 


716( 


-6) 


6.103 


;-6) 


0.9616 


0.9662 


-9 


151 


-4) 





8936 




























0.9802 


0.9832 


-3 


373 


-4) 





9453 


6.024 
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